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INTRODUCTION 


An ancient problem which has intrigued mankind concerns the shape 
of the planet on which men reside. How big is the Earth and what shape 
does it have? To answer this question people developed the science of 
geodesy, the study of the shape and size of the Earth. Traditionally 
this field has often been isolated from other branches of Earth science, 
its goal being narrowly defined as the precise determination of the 
figure of the Earth and distances thereon. As the general shape has 
been known for quite some time, this principal question has been repea- 
tedly answered with successively better accuracy. At this time however, 
geodetic research is undergoing a transition towards greater interaction 
with other Earth sciences such as geophysics and physical oceanography. 

This change is due to new advances in technology, specifically to the 
development of radar altimeters mounted on board Earth orbiting satel- 
lites. The?‘e instruments are capable of mapping the shape of the ocean 
surface and providing accurate estimates of the geoid height. Such an 
advance promises to yield new and exciting information about the Earth. 

As always though in any new field of research there are basic problems 
and questions prior to any new discoveries. In this thesis such fun- 
damental questions are first addressed, and then this information on the 
shape of the ocean surface is utilized to make new inferences about the 
internal structure of the Earth. 

For centuries mathematical geodesists have been able to calculate 
the shape of the Earth. If the internal density distribution is known 
then an appropriate differential equation can be solved for the ellip- 
ticity of the Earth. Alternatively after measurement of the gravity 
field on the surface of the Earth, Stokes' formula can be utilized to 
calculate the geoid. The geoid is defined as a gravitational equipotential 


surface which most nearly coincides with mean sea level. Quite often 
it is termed the shape of the Earth. 

Rather than calculation of this shape, an alternate approach is the 
direct measurement of the shape of the Earth. To a good approximation 
this is possible when utilizing radar altimeters mounted on board 
orbiting space craft. Such an instrument determines the height of the 
satellite above the instantaneous ocean surface. After combination with 
accurate orbital tracking of the satellite motion the height of the 
ocean surface relative to the Earth's ellipsoid is determined. Tech- 
nically the instantaneous sea surface height departs from the geoid by 
up to approximately 1 meter in the deep sea due to oceanographic effects 
such as tides and currents. Since the geoid signal is much larger than 
this, up to 100 meters, these oceanographic influences can be termed 
noise. T+ius satellite altimeters determine the height of the marine 
geoid in the presence of oceanographic noise. For most geophysical 
purposes this is a suitable approximation. 

In April, 1975 an Earth orbiting satellite was launched by NASA. It 
was called GEOS-3, and its mission dedicated to the mapping of the shape 
of the ocean surface. This satellite has provided a wealth of new 
information about the Earth, specifically more than 20 million kilometers 
of measurements of the marine geoid. In order to use this information 
to study the interior of the Earth there are basic questions to be 
answered. Does this Instrument really determine the marine geoid, and 
how well does it function? To partially answer this question it is 
possible to calculate the marine geoid from ship gravity data and make 
comparisons with the satellite altimeter data. This is accomplished in 
chapters 1 and 2. Given measurements of the geoid, what are the 
techniques necessary to interpret the data and make inferences about 
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Internal densities which cause the anomalies? In chapter 3 such methods of 
interpretation are developed. What are the basic results of satellite 
altimetry, what types of geoid anomalies have been detected? A basic survey 
of the various types of geoid anomalies is presented in chapter 4. In 
utilizing geoid data to study the Earth, it is best applied in combination 
with marine gravity data. But in utilizing gravity measurements at sea, 
what corrections are necessary Kecause they were made on the undulating 
geoid? The need for geoid corrections in marine gravity is demonstrated 
in chapter 5. Lastly the question is asked, what new Information about the 
Earth has been learned from the GEOS-3 satellite. In the final chapter a 
study of geoid anomalies over deep sea trenches demonstrates the power of 
marine geodesy and marine gravity in studying the interior of the Earth. 

In order to determine how well the GEOS-3 radar altimeter maps the 
geoid, it is possible to make comparisons with an independently determined 
geoid. In the first chapter this procedure is described for the region of 
the Indian ocean. To calculate the geoid height, a technique using 
differences in gravity and geoid values is utilized. Initially average values 
of free air gravity are obtained for 1® x 1® squares in the region. Then 
values of gravity calculated from the GEM-6 gravitational field model are 
subtracted from the average values to obtain difference gravity data. 
Integration of this difference gravity using Stokes' formula results in the 
calculation of the difference geoid; this is added to the GEM-6 theoretical 
geoid heights to obtain the total Indian Ocean 1® x 1® gravimetric geoid. 

This regional geoid has features such as the Indian Ocean low of -130 m as 
well as a regional high in the Southwest Indian Ocean. 

Utilizing the Indian Ocean geoid as well as gravimetric geoids from 
the North Atlantic and Northwest Pacific Ocean extensive comparisons have 


3 


been made with GEOS-3 radar altimeter estimates of geoid height, this is 
described in chapter 2. A most obvious conclusion from such comparisons 
is the constant offsets between the two data sets. These errors may be 
due to scale errors in the reference ellipsoids, lack of atmospheric 
corrections for the gravimetric geoid, errors in calculation of the 
orbit in the radial direction, or altimeter bias errors. Another observation 
is that altimeter data collected with the intensive mode instrument is less 
noisy than global mode measurements. Even after elimination of long wavelength 
errors, there are still discrepancies between the ^o data sets. In some 
cases this is due to the averaging procedure of the gravimetric geoid con- 
struction. It appears though that after correcting for errors in orbital 
determination and other long wavelength errors, the satellite radar altimeter 
can be utilized for a global mapping of the marine geoid, with resolutions 
greater than are easily achieved with gravimetric geoids. 

After establishing the validity of the radar altinwter measurements, 
it is necessary to develop new techniques to interpret marine geoid data. 

This is accomplished in chapter 3. To demonstrate various properties, ideal 
bodies of constant density are first utilized to calculate the geoid anomaly. 
While such an approach has a certain initial utility, more realistic direct 
methods require the computation of the geoid anomaly due to arbitrary two 
and three dimensional bodies. In developing the formula for a two dimensional 
polygon, the surface integral expression for geoid height is changed to a 
line integral formula via Stokes' theorem and then integrated. For a three 
dimensional polyhedron the volume integral expression for geoid height is 
first converted to a surface integral by using the divergence theorem and 
then to a line integral with Stokes' theorem. To better understand the 
geoid anomaly, the relationship in the wavenumber domain between gravity and 
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geold Is developed. This Indicates that the geold has the same phase 
spectrum as the gravity anomaly but the amplitude spectrum Is amplified at 
long wavelengths. 

A basic need In the study of marine geodesy Is to characterize the 
types of geold anomalies which occur over geological structures. Such 
results In satellite altimetry are Illustrated In chapter 4. Over continental 
margins the geold Increases by a 5 or 6 meter step at the shelf break. At 
deep sea trenches there Is a 10 to 20 meter low which Is 100 to 200 kilometers 
In width. Over the Bismarck sea region there Is a 7 meter high geoid anomaly 
500 km In width. At the Southwest Indian ridge there Is a broad geold high 
at least 500 to 700 km wide and 6 meters In amplitude. Above the Caroline 
Islands the ocean surface has a 5 meter positive anomaly, while at the 
Romanche fracture zone the geold has a 5 meter step similar In shape to the 
bathymetry. 

When measurements of gravity are made at sea they obviously are made 
on the geold. However due to undulations of the geold the free air gravity 
anomaly Is not directly proportional to the gravitational attraction of the 
mass anomalies. If however a correction Is first made for the varying 
distance from the center of the Earth and for the gravitational effect of 
mass between the geold and ellipsoid, then the gravitational attraction of 
the mass anomalies Is obtained. This geold correction becomes necessary when 
studying the medium and long wavelength gravitational ancmalles associated 
with the densities Inside the mantle. This process for correcting the free air 
anomaly 1s discussed In chapter 5. 

Finally In the last chapter geold anomalies over deep sea trenches 
are discussed. Over all active subductlon regions In the Pacific and Indian 
Ocean the geold has a chi!i"acter1st1c shape, thousands of kilometers 1n width. 
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Over the Aleutian Trench beginning about 2000 kilometers seaward of the 
trench axis, the geold has a constant Increase In height as the axis 
Is approached. This slope Is approximately a 20 meter rise over a 
distance of 1800 kilometers. Over the trench Itself there Is an approx- 
imate 13 meter decrease to form a 150 kilometer wide mlnlmun. Iirnied- 
lately landward of the axis there Is an 18 meter Increase until a 
maximum Is reached over the Island arc. Behind the arc the geold 
height gradually decreases by about 10 meters over a 1000 kilometer 
Interval. In the vicinity of all major deep-sea trenches certain features 
form the characteristic geold anomaly; a gradual long Increase as the 
trench Is approached from the ocean basin, a narrow low over the trench 
axis and a geold high but with variable slope landward of the trench. 

For the Tonga- Kermadec region a modeling study was undertaken to 
determine the density structure necessary to produce calculated values 
similar to the observed geold and gravity anomaly. From this work several 
conclusions are possible. Shallower depths landward of the trench repre- 
sent a mass excess which must be mostly compensated. Compensation must 
occur by the presence of a mass deficit somewhere In the underlying mantle. 
In our model a thinner uppermost mantle body provides the compensating 
mechanism. This separation of the mass excess and the underlying mass 
deficit Is the primary cause of the char cterl Stic geold anomaly. Presence 
of the descending lithosphere Is only of secondary Importance In explaining 
the geold anomaly. For the Kermadec Trench the density contrast of the 
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slab with the surrounding mantle Is approximately .02 gm/cm . Further 
studies are necessary for the other trench regions to determine if 
similar principles can be used to model the characteristic geold anomaly. 
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Chapter I 


Detailed 1** x I** Gravimetric Ocean Geold 
and 

Comparison with GEOS-3 Radar Altimeter Geold Profiles 
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Introduction 


A comprehensive set of data on the variation of the sea surface 
height has been obtained In the past three years by means of the GE0S>3 
satellite carrying a radar altimeter. With knowledge of the satellite’s 
position and after corrections for oceanographic effects, one obtains 
along-track variations of the geoldal height over the world's oceans. 

These measured undulations of the geold are extremely valuable for 
geophysical Interpretations because they provide, In part. Information 
about density Inhomogeneities In the earth's upper mantle. The distribution 
of the density Inhomogeneltles reflected in the geold will hopefully 
provide Information on convection within the mantle or on other possible 
driving forces of the earth's lithospheric plates. 

Consequently, it Is of special Interest to determine and study the 
geold on a worldwide scale. In the past, the only way to compute the 
geold over oceanic areas was to apply Stokes' integral to the gravity 
anomalies observed over the earth's surface. Gravimetric geoids over 
oceanic areas have been computed in this way in the western North Atlantic 
(TALWANI, POPPE and RABINOWITZ, 1972), In the Indian Ocean (KAHLE and 
TALWANI, 1973) and In the Northwest Pacific (WATTS and LEEDS, 1977). 

The knowledge of the geoid derived from surface gravity Is essential 
In order to calibrate the radar altimeter, test its performance and 
possibly suggest optimum locations at which the altimeter measurements 
should be performed. Eventually, it will be desirable to compute a 
gravimetric geold independent of satellite altimetry In order to obtain 
the dynamic sea surface height, the sligfit discrepai ''; between the geold 
and sea surface being important in ocean dynamics. 
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In turn, thtre are techniques to recover gravity anomalies from 
6E0S-3 measurements of the sea surface height (e.g. KAHN, SIREY, BROWN 
and A6RAWAL, 1976). The method for obtaining gravity anomalies on the 
basis of satellite radar altimetry can be very useful for determining 
the gravity field In remote ereas such as In the southern hemisphere 
where ship measurements are scarce. In order to test this method and 
define Its accuracy It Is necessary to compare the results with Independently 
obtained gravity values from surface ship gravity measurements In surveyed 
regions. 

It Is the purpose of this paper to compile a new set of T x 1” 
mean free-air gravity anomalies In the Indian Ocean by using all the 
gravity data available up to 1976. compute the corresponding gravimetric 
geold and compare the results with GEOS-3 altimeter derived geold profiles 
(In the Southwest and Northeast Indian Ocean). 

m • 

Sea Gravity Data In the Indian Ocean 
Fig. 1 shows the location of the submarine pendulum observations 
(VENING MEINESZ. 1948; GIROLER and HARRISON, 1957; TALWANI, 1962) as 
well as the ship's tracks along which continuous gravity measurements 
were obtained. 

At present, there are about 100,000 gravity data available In the 
Indian Ocean. 

The data used In this study were obtained by 
a) Lamont-Doherty Geological Observatory of Colunrt>1a University 
during VEMA cruises 18, 19, 20, 22, 24, 29 and 33 and during 
CONRAD cruises 8, 9, 11, 12, 14 and 17. 
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Figure 










b) Institute of Geophysics and Planetary Physics of the University 

of California at Los Angeles and Scripps Institution of Oceanography, 
San Diego, during ARGO cruises "Monsoon" and "Luslad" (HELPER, 

CAPUTO and HARRISON, 1963; CAPUTO, MASADA, HELPER and HAGER, 

1964). 

c) Department of Geodesy and Geophysics, Cambridge University, 
during OWEN cruises 611 through 619, 621 through 626, 110, 111 
and during DISCOVERY cruises 671 and 672 (HYDROGRAPHIC DEPARTMENT, 
1963, 1966; WILLIAMS, 1968). 

d) ESSA (now NOAA) during PIONEER and OCEANOGRAPHER Cruises (U.S. 
DEPARTMENT OP COMMERCE, 1969, 1970). 

e) Bundesanstalt fUr Bodenforschung, Hannover, West Germany, during 
METEOR cruise (PLAUMANN, 1965, personal communication). 

f) Woods Hole Oceanographic Institution, U.S. A., during CHAIN 
cruise 100 (BOWIN, 1973, and BOWIN, personal communication). 

The data obtained prior to 1971 are included In the gravity maps 
(TALWANI and KAHLE, 1975) of the International Indian Ocean Expedition 
(IIOE) Atlas of Geology and Geophysics (UDINTSEV, 1975) as well as in 
the previous Indian Ocean geold paper (KAHLE and TALWANI, 1973). During 
recent Lament cruises (1974, 1977) a considerable amount of data of 
greater accuracy have been accumulated which made it desirable to recompile 
the entire set of 1® x 1° mean gravity data. The main areas of substantial 
Improvements In the surface gravity field of the Indian Ocean are the 
Central Indian Basin including the prominent Indian Ocean gravity low as 
well as the South-West Indian Ocean covei’ing the pronounced South-West 
Indian Ocean gravity high. In the West Pacific we incorporated 1° x 1° 
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mean free-a1r anomalies compiled by WATTS and LEEDS (1977). For the 
Indian peninsjl . we used a set of 1° x 1° anomalies by WOOLARD (1970) 
and values for Madagascar were provided by the International Gravity 
Bureau, Paris (CORON, 1972, personal communication). 

The new set of 1® x 1° free-air anomalies for the Indian Ocean is 
listed in Table 1. The gravity values are referred tc the International 
Reference Ellipsoid (f = 1/297). 

The 1° xl° averages were obtained from revised free-air anomaly 
contour maps by dividing each 1° x 1° square into nine smaller squares. 

The values at the centers of these smaller squares were visually inter- 
polated and averaged to give a value for the 1° x 1° square. In areas 
such as the southwest Indian Ocean the 1° x 1° values were determined by 
averaging individual measurements. 

Method of Geoid Computation 

The following method - proposed and used by many authors including 
STRANGE, VINCENT, BERRY and MARSH (1972) and TALWANI, POPPE, and RABINOWITZ 
(1972) - has been adopted in computing the gravimetric geoid for the 
Indian Ocean: 

1) Calculation of the free-air anomalies based primarily on the satellite 
derived gravity field, in this case the GEM-6 potential field model 
which includes harmonics up to n = 16 (LERCH et al., 1974; SMITH et 
al.. 1976). 

2) Subtraction of the GEM-6 free-air anomalies from the new set of 
1° X 1? mean gravity data (: Difference anomalies 6ag). 

3) Computation of the difference geoid AN (Fig. 2) by applying Stokes' 
integral to the difference anomalies 6Ag. 

4) Summation of the GEM-6 geoid (Fig. 3) and the difference 

geoid aN (: Gravimetric geoid N = Ng^^j g + aN) (Fig. 4). 
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TABLE I 


1“ X 1“ averages of free-air gravity, referred to the International 
Ellipsoid. 
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The advantage of this procedure Is that Stokes' Integration (HEISKANEN 
and MORITZ. 1967) 
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S (y) aag da 
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has to be carried out only over the Indian Ocean because outside this 
area the difference anomalies 6ag are assumed to be zero. 

We recognize that we are Ignoring the difference anomalies outside 
the Indian Ocean, which will givj rise to some errors particularly near 
the boundaries of our area of Integration. These errors however are 
less than 1 meter (STRANGE et a1., 1972). In parts of the Indian Ocean, 
where gravity averages could not be obtained by the procedure described 
above - because of scarcity of data - we assigned the difference anomaly 
6ag to be zero. As we shall see In a later section, this gives rise to 
errors In regions of no ship gravity measurements. This procedure Is 
still superior to methods In which gravity values are truncated after a 
certain radius. 

In Equation (l)aN ■ Difference geold, R ■ mean earth's radius, g » 
mean earth's gravity, y > angular distance between the elenent of area 
da (at which 6&g Is given) and the point of calculation, > Stokes' 
function. Stokes' function S(^) Is defined as 


2 

S(ij<) ■ F(i|() 

s1ni(» 


( 2 ) 
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whera F(^) ■ cosy(» ♦ |s1ni|»(1-5 cos4> - 6 sinVj» - 3 cosi(»ln{sin^ + sln*^)) 

Since S M changes strongly near the effect of a square very close 
to the confutation point cannot be calculated by simply using S( 4 *} with 
being the distance to the center of the square. For this reason the 
squares close to the computation points were subdivided Into smaller 
squares In such a way that the percentage error In calculating S(ij') Is 
below U (TALUANI et al.. 1972). 

The GEM-6 Geold, Difference Geold and 1° x 1° Gravimetric Geold 

In subtracting the 6EM«6 free-air anwnalles (n»16) from the 1* x 1* 
surface data and applying Stokes' Integral to the difference anomalies, 
a difference geold Is obtained which reflects mainly the contributions 
of harmonics 16 < n < 180. The order n>180 corresponds to a wavelength 
of 2 degrees which Is approximately represented by the 1® x 1® mean 
free-air anomalies. As such, the difference geold containing wavelengths 
between about 220 km and 2500 km, can be interpreted as a filtered version 
of the total geold. This range Is of special interest for geophysicists 
because mass Inhomogeneities associated with sea-floor spreading and 
lithospheric motions may produce geold ancxnalles of such wavelengths. 

Fig. 2 shows the confuted difference geold. 

The areal average of the difference oravlty anomalies sag Is 
-.54 mgal and that of the difference geold aN Is -.47m. Thus, the 
systematic error In geold height corresponding to this difference is 
negligible and may not be considered further. The general pattern of 
the difference geold can be characterized by the zero metf contour 
line, indicated by the heavy line on Fig. 2. The entire Northwest Indian 
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Ocean Including the major basins such as the Somali and Arabian Basins 
are associated with a negative difference geold reaching -13 m over the 
Arabian Basin. Another large area with negative difference geoldal 
undulation Is the Northeast Indian Ocean covering the southern Central 
Indian iasin, the West Australian Basin and the Indonesian deep sea 
trench with lowest values over the Sunda Trench (-22 m) and Timor 
Trough (-30 m). 

The northern Central Indian Basin, the major part of the Bay of 
Bengal (Ganges Cone) and the entire Southwest Indian Ocean Is characterized 
by a pronounced difference geold high with highest values over the 
Madagascar Ridge (-*-20 m). While the major actively spreading es are 
not well expressed In the difference geold. the triple Junction of the 
three Indian Ocean ridge branches at 25®S. 70®E clearly stands ou" a 
positive feature ('*■8 m). Also aselsmic ridges Including the Madagascar 
Ridge, the Mascarene Plateau, the Ninetyeast Ridge (-^14 m) and the 
Afanasy Nikitin Seamount chain are associated with a positive difference 
geold. 

The total 1** x 1" geold obtained by adding the difference geold to 
the 6EM-6 geold (Fig. 3) Is shown In Fig. 4. While the difference geold 


23 



reflects wavelengths between 220 and 2500 km, the total geold reveals 
Information with all wavelengths greater than 220 km. As might be 
expected, the long wavelength features are still present In the total 
geold. The prominent Indian Ocean geold low (-130 m) as well as the 
geold highs In the Southwest Indian Ocean and over the western Pacific 
are clearly expressed In the 1® x 1° total gravimetric geold. In addition 
the short and Intermediate wavelength features such as the Mozambique 
Ridge and Basin, the Madagascar Ridge and Mascarene Plateau, the Triple 
Junction, the Ninetyeast Ridge and the Sunda Trench are evident by the 
bending of contours around those structures. The minimum south of India 
has shifted towards the northeast and has decreased In amplitude by 8 m. 
The Somali Basin low appears to connect with the Arabian Basin low. 

GEOS-3 Profiles In the Indian Ocean 
Classically the geold height has been determined by the Stokes' 
Integration of gravity values, as demonstrated In the previous section. 

The altimeter In the GEOS-3 satellite now makes It possible to directly 
measure the elevation of the satellite over the ocean surface. When 
combined with precise orbital tracking one obtains the height of the 
ocean surface relative to an earth ellipsoid. The sea surface is not at 
a constant gravitational potential; nongravitational forces can cause 
slight deviations between the geold and sea surface. Dynamic topography 
due to ocean currents can be 100 cm in amplitude (DEFANT, 1941), tidal 
heights In the deep ocean can also be 50 cm In amplitude (SCHWIDERSKI , 
1977). After correction for such deviations, measurements of the geold 
height may be obtained along the sub satellite track. 
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While this procedure for obtaining geoid heights from altimeter 
measurements is valid, it is subject to certain errors. An extensive 
error analysis has been made for 6E0S-3 measurements in the calibration 
area between Florida and Bermuda (Martin and Butler, 1977). This showed 
an average noise level of .72 m in the intensive mode and 1.81 m in the 
global mode for cumulative altitudes every .1 seconds. Additionally, 
bias values of *5.3 m (intensive mode) and -3.55 m (global mode) were 
discovered for altitude determinations by the altimeter. Accuracy of 
the sea surface height measurements is primarily limited by the orbit 
computations. Long wavelength errors in orbital height can be between 
1-2 m rras and 10 m rms depending upon the method of tracking (H.R. 

Stanley, personal communication). A consistent set of altimeter data 
can be constructed by analyzing measurement differences where two satel- 
lite tracks cross. By least square reduction of such crossover errors 
Rummel and Rapp (1977) were able to effectively eliminate any long wave- 
length orbital discrepancy. In this paper we are primarily interested 
in comparing original GEOS-3 data (from NASA Wallops Flight Center) with 
a gravimetric geoid in the Indian Ocean to see the influence of such 
errors. The orbit numbers and dates for the satellite tracks are 
listed in Table 2. 

Operation of the altimeter is conceptually simple. A radar pulse 
is transmitted downward and the return pulse is received after reflection 
from the sea surface. As there is a finite beam width, the instrument 
I measures the sea surface height over a limited area, with a 14.3 km 

radius for global mode and 3.6 km radius for intensive mode. Operation 
of the altimeter is in one of two modes; global mode or intensive mode 
during which measurements are made every .01024 seconds. This is a sampling 
^ rate of approximately every 65.5 meters. During data processing the sea 

surface height measurements are averaged over a finite time. 
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Table 2 


ORBITS AND DATES FOR SATELLITE TRACKS 


Track 

Designation 

Orbit 

Number 

Day 

Year 

Altimeter 

Mode 

Orbital 

Accurac; 

60141 

1194 

184 

1975 

61obal 

A 

60154 

1237 

187 

1975 

Intensive 

A 

60192 

1647 

216 

1975 

Intensive 

J 

60193 

1653 

216 

1975 

61obal 

D 

60195 

1667 

217 

1975 

Intensive 

D 

60212 

2061 

245 

1975 

Intensive 

0 

60219 

2104 

248 

1975 

Intensive 

0 

60223 

2189 

254 

1975 

Intensive 

0 

60224 

2190 

254 

1975 

Intensive 

D 


A 1 0''‘m rms 
D 3 - 10m rms 
J 3m rms 
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This is a low pass filtering operation. In this paper we only use time 
averaged data usually over a 2.048 second interval. Consequently the 
resultant datum is an average measurement of the sea surface height over 
an area of 14.3 by 13.4 km (global mode) or 3.6 by 13.4 km (intensive 
mode). Knowledge of the beam footprint and the altimeter accuracy is 
necessary in understanding the utility and ultimate resolution of GEOS-3 
measurements. 

The geoid measurements made by the radar altimeter are shown in 
Figure 5 and 6. The sub satellite track is plotted on a Mercator projection, 
the geoid height is plotted perpendicular to the track with positive 
values on the north side. An arbitrary constant offset has been subtracted 
from each profile because we are interested only in the relative changes 
of the geoid. Figure 5 is the Indonesian Island arc region and Figure 6 
is the Southwest Indian Ocean. 

Three GEOS-3 tracks, all approximately perpendicular to the Java 
Trench, are illustrated in Figure 5. Each profile shows a steep long 
wavelength increase in the geoid height of 90 m over a distance of 
approximately 2400 km, with highest values towards the northeast. Not 
being symmetric about the trench axis, the geoid continues to increase 
across the trench until leveling off in the Philippine Sea (Fig. 4). 

Directly over the Java Trench the profiles indicate a geoid low of -10 
m amplitude and 250 km wavelength. Analogous with the gravity low over 
deep sea trenches (WATTS and TALWANI, 1974), this geoid low is primarily 
caused by the bottom topography. 

While the principal features of the Northeast Indian Ocean are deep 
sea trenches and an island arc, the Southwest Indian Ocean is the location 
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Figure 6 

GEOS-3 profiles in the Southwest Indian Ocean. Sub satellite track is 
dotted line, track designation (e.g. G0154) is at start of pass. Geoid 
height is plotted perpendicular to track, positive values are northwards 

and the scale is indicated. An arbitrary constant was removed from each 
track. 28 






Figure 5 

6E0S-3 profiles in the Indonesian Island Arc region. Sub satellite 
track is dotted line, track designation (e.g. G0193) is at start of the 
pass. Geoid height is plotted perpendicular to track, positive values 
are northwards and the scale is indicated. An arbitrary constant was 
removed from each track. 

REPRODUCrrilM ry OF THfi 
ORIGINAL I’A{Ly IS POOR 

29 





of a seismically active spreading center - the Southwest Indian Ridge. 
Topographically high. It Is associated with a 6 m positive geold anomaly, 

500 km wide (Figure 6 and 7, especially tracks G0219 and G0224). Add> 
Itlonally a longer wavelength geold gradient exists. Increasing towards 
the southeast and leveling off at the ridge axis (Figure 4 and 6). 

Comparison of gravimetric geold with GEOS-3 altimeter measurements 

Being a new scientific Instrument, the accuracy and precision of 
the GEOS-3 altimeter must be determined. Obviously a useful test would 
be the comparison of the gravimetric Indian Ocean geold with geold 
heights measured directly by the GEOS-3 radar altimeter. In making the 
comparison characteristics of each method must be understood. The 
gravimetric geold contains Information from the GEM-6 gravity model up 
to n = 16, the higher harmonics (to n = 180) being dependent upon the 
accuracy and availability of ship gravity data. Resolution of features 
with wavelengths less than about 200 km cannot be expected In the 
gravimetric geold due to our averaging procedure over 1° x 1® squares. 

With the GEOS-3 altimeter sea surface features with wavelengths of 1.31 
km can technically be resolved (with average values every 655 m), however 
with the averaged data we utilized resolution Is limited to 28 km. Also 
Inaccurate orbit determinations result In constant and very long wavelength 
errors In the altimeter measurements. Consequently we only make comparisons 
of geold Information with wavelengths shorter than several thousand km. 

GEOS-3 measurements of geold height are referred to an ellipsoid 
with flattening f = 1/298.255 and semi-major axis a = 6378145 m, thus 
the gravimetric geold was converted to this ellipsoid prior to com- 
parison. Linear Interpolation of the gravimetric geold with the nearest three 
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geo1d values was used In order to compute the gravimetric geold height 
along the sub satellite track. Results of this comparison are In Figure 
7, each track location Is plotted In Figure 5 or 6 and the track designation 
(e.g. G0195) Is at the start of the pass (time ■ 0 seconds). 

Obviously the radar altimeter cannot measure the geold over land, 
thus such measurements have been deleted (e.g. Madagascar). GEOS-3 
measurements In Figure 7 are Indicated by the thin line and the gravimetric 
geold by the thick line. An obvious feature of this comparison Is the 
offset values of up to 25 m, probably due to error In orbit determination. 
These profiles Indicate a great similarity (except for a constant shift) 
between GEOS-3 data and the gravimetric geold wherever good gravimetric 
control exists (see e.g. G0242). In some areas (south of 39**S In SW 
Indian Ocean) only few sea gravity measurements exist. Consequently, In 
this area the gravimetric geold contains little Information for wavelengths 
shorter than n > 16 (2500 km). This Is seen In tracks G0224 and G0219 
between 0 and 120 seconds over the Southwest Indian Ridge. The GEOS-3 
data Indicate a clear anomaly over the ridge while the gravimetric geold 
does not. As there Is poor gravity control In this area and the anomaly 
Is seen on both tracks, we conclude the gravimetric geold Is Incorrect 
along this part of profiles G0224 and G0219, whereas GEOS-3 measures the 
true geold anomaly over the Southwest Indian Ridge. 

In finer detail, the GEOS-3 altimeter has greater resolution than 
the 1® X 1® gravimetric geold (of this paper), due to a smaller sampling 
Interval for the altimeter. Over the Java Trench (tracks G0195, G0193 
In Figure 7) the effect of the averaging procedure for the calculation 
of the gravimetric geold Is clear: GEOS-3 data indicate a 250 km wide 
low, 10 m deep. The gravimetric geold barely defines the low. Due to 
the Information content In wavelengths shorter than about 200 km, it is 
understandable that the GEOS-3 altimeter with a smaller sampling interval 

would measure the geold more accurately In this region. 

31 



COMPARISON OF GRAVIMETRIC GEOID WITH GEOS-3 ALTIMETER MEASUREMENTS 







Figure 7 


Conclusions 

On the basis of marine gravity measurements In the past three 
decades, we have compiled and listed In Table 1 1° x 1” mean free-air 
gravity anomalies for the Indian Ocean. These values are useful In 
geold computations and as a test for techniques of gravity recovery from 
GEOS-3 geold measurements. Utilizing these r x r averages we have 
computed a difference and total gravimetric geold for the region. 

The difference gravimetric geold contains Information with wavelengths 
between about 200 km and 2500 km, such wavelengths being determined by 
lateral density Inhomogeneities within the crust and upper mantle. 

Anomalies of -18 m In the difference geold exist over the Arabian Basin, 
another large area with negative values Is the Northeast Indian Ocean 
with lowest values over Sunda Trench (-22 m) and Timor Trough (-30 m). 

Over the Madagascar Ridge there Is a 20 m difference geold high, another 
positive feature Is the triple Junction of the three Indian Ocean ridges, 

■•■8 m amplitude. Aselsmic ridges such as Madagascar Ridge, Ninetyeast Ridge 
and Afanasy Nikitin seamount chain are also associated with a positive 
difference geold. 

The total gravimetric geold contains Information on all wavelengths 
down to about 200 km; thus It has both long and Intermediate wavelength 
features. Over the Central Indian basin Is the most prominent long 
wavelength feature, the -130 m Indian Ocean geold low. Steep geoldal 
gradients exist over the Indonesian Island Arc up to a flat region of 
the geold In the Phllllpine Sea. Another level portion of the geold 
exists south of the Crozet Plateau. Intermediate wavelength anomalies 
In the total geold are over the same features as In the difference geold. 
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6E0S-3 profiles of geold height show anomalies across topographic 
f features In the Indian Ocean. Over the Southwest Indian Ridge there Is 

a 6 m ? 'itive geold anomaly, 500 km wide. Across the Java Trench there 
Is a steep Increase In the geold of 90 m over a distance of 2400 km, 
with highest values towards the northeast. Directly over the Java Trench 
there Is a -10 m geold low, 250 km wide. If this relative low Is removed, 
an overall geold high remains. This Is an Interesting result because It 
might shed some light on the density Inhomogeneities associated with the 
descending Indian lithospheric plate. In part, this pronounced geold 
high - verified In the 6E0S-3 profiles - can be explained by the positive 
density contrast of the cold lithosphere with respect to the less dense 
adjacent asthenosphere. On the basis of gravity data Watts and Talwani 
(1974) concluded this effect Is not the only component because unrealistically 
high density values would have to be postulated to explain both the 
"outer gravity high" and the overall high. It Is further Interesting to 
note that In our 6E0S-3 profiles an "isolated" outer geold high does not 
exist. The Increase In geold height seaward of the Indonesian Trench 
appears to be a portion of the overall long-wavelength high - only 
Interrupted by the relative small -wavelength low over the trench proper. 
Another possible source for part of this long-wavelength geold high may 
be associated with the downwarped Isotherms caused by the cold sinking 
lithosphere. If the temperature Is lowered at the Olivine-Spinel transition 
zone (at depths between 300 and 400 km) this phase change migrates 

3 

upwards, thus providing a positive density contrast of about 0.2 gm/cm 
(BOTT, 1971; RINGWOOO, 1976). To explain this long wavelength geold 
anomaly at the Java Trench, a detailed modeling study should be undertaken; 
only tnen will the cause of this anomaly be better defined. 
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Co(np«r1sons of GEOS-3 data with tha gravimetric geold show a great 
similarity where good gravimetric control exists. Due to the shorter 
sampling Interval (14 km versus about 100 ion averaging for the gravimetric 
geold) for the altimeter the GEOS-3 satellite can better detect short 
wavelength features such as the low over the Java Trench. On the basis 
of these conparlsons we conclude the GEOS-3 altimeter Is an accurate and 
highly useful Instrument for mapping the geold at sea. 




CH APTER 1 ACKNOWLEDGEMENTS ; 

Principal authors for Chapter 1 were: 1) Hans-Gert Kahle of 
Instltut fur Geophysik, ETH-Honggerberg, CH-6093 Zurich, Switzerland 
and Lamont-Doherty Geological Observatory of Co1unA>1a University, 
Palisades, New York 10964, U.S.A.; 2) Michael Chapoian and Hanik Talwani 
of Lamont-Doherty Geological Observatory of Coluir^la University, 
Palisades, New York 10964. 

We are grateful to A. B. Watts and J. R. Cochran for critically 
reading the manuscript and for helpful suggestions for Its In^rovement. 
This research was partially sponsored by Office of Naval Research, 
Contract N0OQ14-75-0210; La>nont-Ooherty Geological Observatory of 
Columbia University, Contribution No. 2672 and Instltut fur Geophysik, 
ETH Zurich, Contribution No. 176. 


36 




REFERENCES 



8ott» 1971. Tht mantis transition zona as a possible source of 

global gravity anomalies, Earth and Planetary Science Letters, 11, 

28. 

BoMln, C., 1973. Origin of the N1noty*East Ridge from studies near the 
equator, J. Geophys. Res., 78, 6029. 

Caputo, N., Hasada, R., Heifer, M.D., and Hager, C.L., 1964. Gravity 
measurements In the Atlantic, Pacific and Indian Oceans, May 1962 - 
August 1963, R/y ARGO. Interim Report, Univ. of Calif, Inst, of 
Geophysics and Planetary Physics. 

Defant, A., 1941. Quantitative Untersuchungen zur Statik und Dynamlk 
des Atlantlschen Ozeans, Ole relative Topographic cinzelner 
Druckf lichen In Atlantlschen Ozean, Deutsche Atlantlsche Exped. 

Meteor, 1925-1927, Miss. Erg.. VI, 183. 

Girdler, R.W. and Harrison, J.C., 1957. Submarine gravity measurements 
In the Atlantic Ocean, Indian Ocean, Red Sea and Mediterranean Sea. 
Proc. Roy. Soc., A 239, 202. 

Helskanen, W.A. and Moritz, H., 1967. Physical Geodesy, W.H. Freeman 
and Co., San Francisco, Calif. 

Heifer, M.O., Caputo, M. and Harrison, J.C., 1963. Gravity measurements 
In the Pacific and Indian Oceans, MONSOON Expedition 1960-196i. 

Interim Rept., Univ. of Calif., Inst, of Geophysics and Planetary 
Physics. 

Hydrographic Department, Admiralty, 1963, 1966: Bathymetric, magnetic 

ar.d gravity Investigations H.M.S. OWEN 1961-1962, 1962-1963, Admiralty 

Marine Scl. Publ., no. 4, parts 1, 2; no. 9, parts 1, 2. 

Kahle, H.G. and Talwani, M. 1973. Gravimetric Indian Ocean Geold. Z. 

fflr Geophysik, 39, 167 and 491. 

37 


Kahn, M.O. . Sirey, J.W., Brown, R.D. and Agrawal, B.B., 1976. Gravity 
anomaly determination by means of GEOS-3 altimetry data. (Abstract 
only) EOS, 57, 901. 

Lerch, F., Wagner, C.A., Richardson, F.A., Brown, F.E., 1974. Goddard 
Earth Morf^ls 5 and 6 Rep. X-921-74-145, Goddard Space Flight Center, 
Greenbelt, Md. 

Martin, C.F. and Butler, M>L., 1977. Calibration results for the 
GEOS-3 altimeter, NASA Contractor Report CR-141430. 

Ringwood, A.E., 1976. n se transformations In decending plates and 
Implications for mantle dynamics, Tectonophysics, 32, 129. 

Rummel, R. and Rapp, R., 1977. Undulation and anomaly estimation 

using GEOS-3 altimeter data without precise satellite orbits. Bull. 
Geod. , 51, 73. 

Schwidsrski , E.W., 1977. Ocean tides and GEOS-3, Trans., A.G.U., 58, 

370 (Abstract only). 

Smith, D.E., Lerch, F.J., Marsh, J.G., Wagner, C.A., Kolenklewicz, R., 
and Kahn, M.A., 1976. Contributions to the National Geodetic 
Satellite Program by Goddard Space Flight Center, J. Geophys. 

Res. , 81 , 1006. 

Strange, W.E., Vincent, S.F., Berry, R.H. and Marsh, J.G., 1972. Detailed 
gravimetric geold for the United States, The Use of Artificial 
Satellites for Geodesy, S. Henri ksen, A. Mancinl, and B. Chovltz 
(eds.), Washington, O.C., 169. 

Talwani, M., 1962. Gravity measurements on H.M.S. ACHERON In South 
Atlantic and Indian Oceans, Geological Society of America Bull., 

73, 1171. 

Talwani, M., Poppe, H. and Rablnowitz, P., 1972. Gravimetrical ly 
determined geold In the western North Atlantic. Iji: Sea surface 

topography from space, NOAA Tech. Rept. ERL-228 - AO ML 7-2, 34p. 

38 


Talwani, M. and Kahle, H.-G., 1975. Maps of free-air gravity anwnalies 
in the Indian Ocean, In: The International Indian Ocean Expedition 
Geological - Geophysical Atlas of the Indian Ocean, G. Udintsev 
(ed.), Moscow, 87. 

Udintsev, G. (ed.), 1975. The International Indian Ocean Expedition 
Geological - Geophysical Atlas of the Indian Ocean, Moscow. 

U.S. Department of Commerce, Environmental Science Services Administration, 

1969. International Indian Ocean Expedition USC & GS Ship PIONEER, 
1964, V. 1 & 3. Washington. 

U.S. Department of Commerce, Environmental Science Services Administration, 

1970. Technical Report ERL 152-00 4. Global Expedition USC & GS 
Ship OCEANOGRAPHER - 1967 Cruise Narrative and Scientific Program. 
Boulder. 

Vening Meinesz, F.A., 1948. Gravity expedition at sea 1923-1938, v.IV, 
Publ. Neth. Geod. Commission, Delft. 

Watts, A. 8. and Talwani, M., 1974. Gravity anomalies seaward of deep- 
sea trenches and their tectonic implications, Geophys. J. Roy. 

Astr. Soc., 36, 57. 

Watts, A.B. and Leeds, A., 1977. Gravimetric geoid in the Northwest 
Pacific Ocean, Geophys. J. Roy. Astr. Soc., 50, 249. 

Williams, C., 1968. R.R.S. DISCOVERY cruise 16 Indian Ocean 1967, 
bathymetric, gravity and magnetic data report. 

Woollard, G.P., 1970. 1“ x 1" square free-air gravity averages in 
India, Publ. Hawaii Inst, of Geophys. 


39 


Chapter II 


Comparison of Gravimetric Geoids with GEOS-3 Altimetric Geoid 
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INTRODUCTION 

With the aid of a radar altimeter mounted on board the GEOS-3 
satellite It Is now possible to rapidly determine the shape of the ocean 
surface. This Information can be used to provide an estimate of the 
marine geold. However, because the altimeter Is a new Instrianent, It Is 
of considerable Interest to compare such measurements with other types 
of data. Thus we examine how well 6E0S-3 estimates of geold height 
compare with data from Independently determined gravimetric geolds. 

GRAVIMETRIC GEOIDS 

On land the geold may be constructed using astrogeodetic n«thods; 
or alternatively utilizing gravimetric methods and Stokes' Integration 
of the measured values of gravity. At sea direct determination of 
deflection of the vertical is a difficult task, consequently only about 
20 measurements have ever been made (Von Arx, 1966). Due to this dif- 
ficulty, only with measurements of gravity at sea can geoid computations 
be performed. This technique for oceanic geoid construction was demon- 
strated in the western North Atlantic Ocean (Talwani, et al, 1972). 
Subsequent studies have involved construction of oceanic geolds in the 
Indian Ocean (Kahle and Talwani, 1973), Northwest Pacific (Watts and 
Leeds, 1977) and an updated Indian Ocean geoid (Kahle, Chapman and 
Talwani, 1978). Worldwide geoids have been constructed by Heiskanen 
(1957) and Marsh and Vincent (1974). 

Basically our technique of geoid construction involved a series of 
computational steps: 
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1. Compilation of marine gravity data and averaging over a specified 
area - in this case 1“ x 1® averages. 

2. Calculation of free air anomalies on the basis of a certain 
gravitational potential model; our computations utilized the 
GEM-6 model (Smith et al., 1976). 

3. Subtraction of the calculated free-air anomalies (GEM-6) from 
the areal averages of measured gravity; this is the set of 
difference gravity anomalies. 

4. Application of Stokes' integral to these difference gravity 
values to obtain the difference geoid. 

5. Addition of this difference geoid to the geoid of the gravi- 
tational potential model (GEM-6) to obtain the total gravi- 
metric geoid. 

Computationally this procedure is efficient because it eliminates 
the need for calculation of Stokes' integral over the entire earth. As 
the integration is done only over the area of study, this is equivalent 
to utilizing GEM-6 gravity values outside the area of computation. Some 
errors can be introduced, especially near the boundaries; however, 
such errors are less than 1 meter (Strange et al., 1972). 

Due to this procedure there are several inherent characteristics of 
these gravimetric geoids. Outside a certain region the values of 
gravity due to a satellite model are utilized, consequently the long 
wavelength (roughly N<^16) components in these geoids are determined by 
that particular gravity model. In our case the long wavelength com- 
ponents of the gravimetric geoids are determined by the GEM-6 field. 
Another property of these geoids results from the initial procedure of 
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averaging gravity values over a finite region. Averaging over a certain 
interval and subsequent decimation of data corresponds to a low pass 
filter operation with a gradual cutoff band. In our gravimetric geoids 
this procedure would eliminate most wavelengths shorter than 2** which is 
twice the sampling period. At a latitude of 23° then one could state 
that our 1° x 1° gravimetric geoids do not contain much information in 
the wavelengths shorter than about 200 km. Another characteristic of 
the present marine geoids is that in some regions where there are few 
ship tracks there is difficulty in obtaining average gravity values. As 
Stokes' function is most sensitive to adjacent locations, it would be 
expected that in regions where gravity averages are poorly determined, 
the geoid would likewise not be accurately calculated. 

In constructing these geoids there were no corrections made for the 
effect of mass of the atmosphere external to the geoid (Moritz, 1974). 

In ignoring this problem, errors of several meters in geoid height are 
possible (Rapp, 1975). Such errors would be manifest as very long 
wavelength differences or constant offsets in geoid height. 

Prior to our comparison study, the gravimetric geoids were first 
transformed to be relative to a best fitting ellipsoid with flattening 
1/298.255. During the transformation only the J 2 and terms were 
altered; there was no change in the zero'th order undulation. Theore- 
tically then the radius of this new reference ellipsoid corresponds to 
the actual radius of the earth. GEOS-3 data is referred to an ellipsoid 
with flattening 1/298.255 and radius 6378142 m. There is a possible 
scale difference between the reference ellipsoids for gravimetric geoids 
and altimetric geoids. They do, however, have identical values of 
flattening. 
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For purposes of comparison in this study we utilize gravimetric 
geoids in the Indian Ocean (Kahle, Chapman and Talwani, 1978) as shown 
in Figure 1, Northwest Pacific (Watts and Leeds, 1977) in Figure 2, and 
North Atlantic (Talwani and Leeds, in preparation) in Figure 3. 

GEOS-3 ALTIMETER 

Several characteristics of the 6E0S-3 radar altimeter have a 
bearing on our study. On board the satellite there are two transmit- 
ters, one designated the global mode and the other termed intensive 
mode. Both measure the height of the satellite above the sea surface 
every .01 seconds. In practice it has been found that the global mode 
results in a signal with higher noise; consequently this mode was little 
used after the initial stages of the mission (H.R. Stanley, personal 
communication). During computer processing of this data an average 
measurement is determined. For the low rate telemetry format this is a 
2.048 second average and for the high rate telemetry format a 3.277 
second average. At an average ground track speed of 6.55 km/sec this 
results in a measurement every 14.3 km (low rate) or 21.5 km (high 
rate). Due to this averaging procedure most wavelengths shorter than 
28.6 km (low rate) or 43 km (high rate) will be eliminated from our 
GEOS-3 measurements. 

After the radar altimeter measures the altitude of the spacecraft, 
orbit calculations are utilized to locate the satellite relative to the 
center of mass of the earth. During the orbit computation a satellite 
derived gravity field is used; any errors in the assumed geopotential 
coefficients will generate uncertainty in the radial position of the 
satellite. In a comparison study of the SKYLAB radar altimeter and 
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Figure 




Figure 2 










satellite geoids, an error analysis of orbital uncertainty was performed 
by Marsh et al. (1976). Due to possible errors In the quantity GM there 
was a mean radial uncertainty of 2.8 m In the SKYLAB orbit. Possible 
errors In geopotentlal coefficients and station coordinates generated 
radial orbit uncertainty of approximately +1 m about the n»an value. 

This uncertainty In radial position of SKYLAB was slowly varying over a 
given orbit. Certainly for the 6E0S-3 satellite an orbital error ana- 
lysis would Indicate different numerical values; however, similar Infer- 
ences can still be made. Any errors In the quantity GM, coefficients of 
the geopotentlal, and station coordinates will all result In mean radial 
uncertainties and slowly varying errors In the radial position of GEOS-3. 
What causes additional confusion with this satellite (GEOS-3) Is that 
several different gravitational field models were utilized for the orbit 
computations; for the data In this paper It was a GEM-8 (Wagner et al., 
1976) or an NWL model. Another possible source of uncertainty Is caused 
by any bias In altimeter range data. This Is discussed In detail by 
Martin and Butler (1977). 

After altitude measurement and orbit calculation, the sea surface 
height Is determined. To obtain an estimate of geold height, tidal 
corrections should be made. Our GEOS-3 data were corrected with the 
tide model of Hendershott (1973). 

COMPARISON OF GRAVIMETRIC GEOIDS WITH GEOS-3 ALTIMETER 

Using GEOS-3 altimeter estimates of geold height, we have made 
comparisons with the gravimetric geoids presented in the previous 
section. The gravimetric geoids were Interpolated at the point of 
altimeter measurement; interpolation was linear relative to the nearest 
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three gravimetric geo1d values. A track chart of the location of 
altimeter measurements, and the extent of the gravimetric geolds Is 
shown In Figure 4; a listing of orbit numbers for each track Is In 
Table 1. 

In Figure 5 there Is a comparison of data with the Indian Ocean 
geold, the altimeter data Is plotted as a function of time of acquisi- 
tion from the first data point. The first obvious fact from this com- 
parison Is that GEOS-3 data often differs from the gravimetric geold by 
a constant level. RMS differences between the two sets of data are In 
the range 3.69 m to 13.27 m as listed In Table 2. 

In analyzing the causes for this discrepancy there are a variety of 
possible causes. These Include bias In altitude determination (Martin 
and Butler, 1977), radial orbital uncertainty, scale errors In the 
geold, and long wavelength errors In either the gravimetric geold or 
radial orbital position. To test for the effect of bias In altitude 
measurement, corrections for this error were made using the results of 
Martin and Butler (1977). After correction for the effects of bias 
utilizing the results of Martin and Butler, the rms differences were 
recomputed and results listed In Table 2. In general this bias cor- 
rection often Increases the discrepancies. This shows that our com- 
parisons of the GEOS-3 values are mostly lower than the gravimetric 
geold, and the bias correction enlarges this difference. Due to small 
uncertainties In the geopotential coefficients, value of GM, and trac- 
king station coordinates, there may be uncertainty In the radial orbital 
position of GEOS-3. Such errors would be manifest as either constant or 
slowly varying radial uncertainty. This has been proven In a quan- 
titative manner for the SKYLAB altimeter by Marsh et al. (1976). Any 
radial error in orbit determination can generate apparent constant 
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Table 1 


Track no. Orbit Telemetry mode Transmitter mode 


INDIAN OCEAN 


G0118 

646 

Low 

Intensive 

G0275 

461 

Low 

Intensive 

G0142 

1200 

Low 

Global 

GDI 82 

1568 

Low 

Intensive 

G0106 

362 

High 

Intensive 

G0112 

584 

Low 

Intensive 



NORTHWEST PACIFIC 


G0242 

2037 

Low 

^1tensive 

G0230 

1795 

Low 

Intensive 

G0235 

416 

Low 

Global 

G0266 

430 

Low 

Global 



NORTH ATUNTIC 


G0329 

3245 

Low 

Intensive 

G0039 

184 

Low 

Global 

G0096 

325 

Low 

Global 

G0069 

398 

High 

Intensive 

G0031 

210 




NORTHWEST PACIFIC 


G0846 

1724 

Low 

Intensive 

G0855 

1616 

Low 

Intensive 

G0565 

2028 

Low 

Intensive 
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Comparison of Indian Ocean gravimetnc geoid with GEOS-3 
estimates of geoid height. Geoid height is plotted as a 
function of time of acquisition from first point. The sub 
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Table 2 


RMS discrepancies between gravimetric geolds and GEOS-3 estimates of geold 


Track no.* RMS difference RMS difference A priori orbital error 

after bias cor- 
rection 


INDIAN OCEAN 


G0106 

3.69 meters 

4.78 meters 

10* meters 

G0112 

6.99 

12.26 

K 

G0118 

13.27 

18.37 


G0142 

6.67 

10.01 

K 

G0275 

10.84 

16.02 

10* 

G0182 

11.15 

6.46 

3 



NORTHWEST PACIFIC 


G0242 

17.87 

15.87 

3 

G0230 

18.00 

23.61 

3t10 

G0235 

5.86 

11.07 

K 

G0266 

6.82 

10.27 

10 



NORTH ATLANTIC 


G0031 

18.99 



G0039 

23.93 

37.23 

K 

G0069 

13.18 

15.65 

10 

G0096 

13.18 

15.65 

3-10 

G0329 

3.07 

8.05 

3 


*0rb1t numbers listed In Table 1 
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offsets between the gravimetric and altimetric geolds. If the actual 
orbit error is slowly varying, then over the short segment during which 
altimeter data is acquired it may appear to have a constant offset or 
tilt which in reality is only one part of a mucn longer wavelength error 
in radial position of the satellite. A complicating factor with GEOS-3 
is that several different sets of geopotential coefficients have been 
utilized for the orbit calculations. These gravity models are either 
GEM-8 (Wagner et a1., 1976), GEM-10 (Lerch et al., 1977), or a NWL 
model. In properly utilizing GEOS-3 data it is mandatory that all 
orbits be calculated with only one geopotential model. This would 
provide a consistent data set. 

During orbit calculations by NASA a priori estimates of rms orbital 
errors were made. These are listed in Table 2. For example, track G0106 
has an a priori radial uncertainty of greater than 10 m rms while the 
rms discrepancy between the gravimetric and altimetric geoid is 4.78 m. 

In other cases the actual rms difference exceeds the prior estimates. 
Additionally constant offset discrepancies between the two data sets may 
occur due to scale errors in the gravimetric geoid or ellipsoid for 
GEOS-3 data and to the lack of atmospheric corrections with the gravi- 
metric geoids. In this type of comparison study it is virtually impos- 
sible to isolate the cause of the observed constant offsets between 
GEOS-3 geoids and gravimetric geoids. They may be caused by any or all 
of the aforementioned sources of error. These large constant offsets do 
indicate the need for high quality orbital computations with a con- 
sistent gravity model, atmospheric corrections in gravimetric geoids, 
elimination of errors in reference ellipsoids, and analysis of internal 
consistency of the data. 
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By analyzing altimeter data at crossover locations and least square 
reduction of crossover errors, Rummel and Rapp (1977) were able to 
obtain rms crossover discrepancies of .78 m for Intensive mode opera- 
tion. Our comparisons Indicate the need for this type of crossover 
analysis and error reduction. 

Comparisons of the altimeter data and the Northwest Pacific gravi- 
metric geold are shown In Figure 6. In this the Influence of the 
transmitter mode Is evident. Tracks 60242 and 60230 are relatively 
smooth and were collected In the Intensive mode, tracks 60235 and 60266 
have high frequency noise and were made In the global mode. In general 
all of the figures show that the Intensive mode has less high frequency 
noise than the global mode of acquisition. 

Another aspect of our comparison study Involves the question of 
what are the shortest wavelengths In the geold and how well does the 
gravimetric geold and altimeter data record such wavelengths. To under- 
stand the short wavelength (less than several hundred kilometers) com- 
ponents of the geold, It Is necessary to examine their origin. Any 
geold anomaly Is due to the anomalous potential caused by a mass hetero- 
geneity. On the basis of potential theory It would be expected that 
short wavelength components of the geold would be caused by the nearest 
mass anomalies; this would be topography of the ocean floor and moho 
which are large, adjacent, mass Inhomogeneities. To see this mathematically 
consider the relationship between the Fourier transform of gravity and 
bathymetry, the admittance. 

Z(k) = 






Cg(k)] //" [b(k)] 


( 1 ) 
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Comparison of Northwest Pacific gravimetric geoid with 
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labeled by the Lamont track number, gravimetric geoid 
heights are the smooth line. 
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For an Airy model of two dimensional crust this function would be 
(McKenzie and Bowin, 1976) 

Z(k) ■ 2.G -yo^) e-''’' (1 - e''“) (2) 

where G Is the gravitational constant 
and^j.,y^, are density of crust and water respectively 
k Is wavenumber In radians/ km 
d Is depth of the water layer 
t Is the thickness of crustal layer 

Utilizing the transfer function between the Fourier transform of gravity 
and geold (Chapman, 1979; A. Leeds, personal communication) we can 
obtain the admittance between the Fourier transform of geold and bathy- 
metry for two dimensional Airy Isostasy: 

■ .7^ (A -^w> (3) 

where 5^N) Is the Fourier transform of 
geold height 

Y 1s normal gravity : 980 cm/sec 

What this equation shows Is that topography of the ocean floor will give 
rise to undulations of the geold which have Identical frequencies but 
amplitude decreasing with an increase In frequency at high wavenumber. Any 
other model of compensation would also have Identical frequencies but a 
different amplitude function. This relationship 1s Important in our 
comparison study because It indicates that topographic features will 
generate undulations of the geoid with similar frequencies. In the case 
of the Hawaiian ridge this has been proven in a quantitative manner by 
Watts (1978). 


57 


As an example of this relation between bathymetry and geoid con- 
sider the Aleutian trench (Figure 7). At the trench axis there is a 
topographic depression which is approximately 100 km wide, it has a 
gravity anomaly of similar width. Assuming two dimensionality it is 
possible to compute what the geoid anomaly is from the gravity anomaly 
(Taiwan! et al, 1972; A. Leeds, personal communication). 

N(x) » (^) g(x') log lx - x'l dx' (4) 

This is a convolution integral and is the spatial equivalent of the 
transfer function between the Fourier transform of gravity and geoid 
(Chapman, 1979). Utilizing this formula and the observed gravity across 
the Aleutian trench, a two dimensional geoid profile has been computed 
and is also shown in Figure 7. Additionally a GEOS-3 profile and 
corresponding gravimetric geoid is shown. From this it is clear that 
above the trench there is a two dimensional geoid minimum with similar 
frequencies as the bathymetry. In both amplitude and wavelength this is 
identical to the observed minimum in GEOS-3 data. However, the trench low 
as seen in the 1® x 1® gravimetric geoid is much broader and shallower. 

On the basis of the relation between geoid and bathymetry having iden- 
tical frequencies and our computation of the two dimensional geoid, it 
is clear that above the narrow Aleutian trench there is a geoid anomaly 
with similar frequency components. As the GEOS-3 altimeter records a 
signal similar in shape to the two dimensional geoid, it is reasonable 
to assume that the altimeter is faithfully measuring the geoid signal. 

In the 1® X 1° gravimetric geoid the trench low is much broader and 
shallower; in comparison to the GEOS-3 data and two dimensional geoid, 
the gravimetric geoid has less high frequencies and more energy in the 
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Geold height and bathymetry across the Aleutian trench. 
Bathymetry profile is at bottom and is from a ship track* 
the gravity profile is also from Conrad cruise 1109. Utilizing 
this gravity profile which is projected normal to the ^ 

trench axis, and formula 4, a two dimensional neoid, was 
shown. This compares well w ih the 
GEOS-3 altimetry data; however the 1® X 1® grav ...etric 
poid has a wider and shallower low above the trench axis. 

0^ the projected ship track and 

GEOS-3 profile is also shown. 
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medium frequency (wavelengths approximately 200 km) range. This is 
aliasing of the geoid signal and is due to our procedure in geoid con- 
struction of averaging gravity values over 1® squares. 

This analysis of geoid height over the Aleutian trench indicates a 
limitation of 1® x 1® gravimetric geoids. In certain regions there can 
be energy in the geoid for wavelengths shorter than are resolved by 
averaging over 1® squares. In such areas it is necessary to construct 
gravimetric geoids by first averaging gravity values over smaller 
regions, perhaps 10' or 5' squares. Prior to averaging these values, it 
is possible to estimate how much high frequency information exists in 
the gecid. This can be done by multiplication of the Fourier transform 
of bathymetry and the admittance function. This gives an estimate of 
what the geoid heights would be at the shorter wavelengths. Another 
technique to estimate the high frequency geoid heights would be to first 
compute the Fourier transform of gravity values. After multiplication 
by the transfer function for a plane earth, this yields the Fourier 
transform of geoid height. This provides an estimate of the geoid 
height at the highest frequencies. 

Data from the North Atlantic region are compared in Figure 8. 

Again the same features are noted as before, constant offsets and noisy 
altimeter data when the global mode is utilized. However, track G0096 
exhibits another feature: there is a tilt of the altimeter measurements 
relative to the gravimetric geoid. Because this is such a long wave- 
length difference and does not appear in the other regions of the North 
Atlantic, we suspect it is due to very long wavelength errors in the 
computed satellite orbit. For this reason, Rurimel and Rapp (1977) in 
correcting GEOS-3 data removed both a long wavelength orbital tilt in 
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addition to constant offsets. Such a procedure Is both justified and 
necessary In order to obtain good geold estimates from the GEOS-3 
altimeter. 

In an effort to overcome such long wavelength differences, four tracks 
of 6E0S-3 altimeter data were adjusted by requiring that they agree with 
each other at crossover locations and with the GEM-6 geold. To do this 
adjustment the technique of Rummel and Rapp (1977) was utlllz this 
minimizes In a least square sense discrepancies due to bias and tilts. 

A compar'ison of these adjusted GEOS-3 geold estimates with the Northwest 
Pacific gravimetric geold Is shown In Figure 9. For each track the 
GEOS-3 data and gravimetric geold are both shown, directly above this 
the difference In height between the two Is shown at a different scale. 

Thus for track G0266 agreement Is quite good except at the Bonin trench; 
the rms discrepancy Is 2.07 m (Table 3). For all of these tracks major 
disagreements of up to 10 m occur above the trench systems such as the 
Bonin, Mariana, Ryukyu, and Philippine trench. As discussed previously 
for the Aleutian trench, this discrepancy Is due to inadequate reso- 
lution of the gravimetric geold. 

Other broad regions of difference occur; on track G0565 (Figure 9) 
southeast of the Mariana trench there Is a broad disagreement of up to 5 
meters. Because this area does not have the same quantity of gravimetric 
data as other regions, we suspect this Is due to errors In the gravi- 
metric geold. When making detailed comparisons though. It becomes 
exceedingly difficult to ascertain whether there are slight errors In 
orbital computation or geold calculation, or simply one could be seeing 
the effect of transient sea surface topography in the altimetry data. 

With gravimetric geolds what Is needed are finer resolution grids and 
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Figure 9 
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Table 3 


RMS discrepancies between gravimetric geolds and adjusted GEOS-3 data 

In Northwest Pacific 


Track no. RMS difference before RMS difference after 
adjustment adjustment 


60846 

12.34 meters 

2.67 meters 

60855 

4.82 

1.69 

60565 

4.rj 

2.36 

60266 

10.27 

2.07 
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detailed error analysis of ti»e1r quality: In orbit determination the 
best calculations should be made to reduce radial errors to less than . 
meter. 

The primary miss un of the 6E0S-3 satellite was to determine Infor- 
mation about the gravitational field of the earth; another question then 
arises as to what wavelength Information will be best determined by 
satellite altimetry. For very long wavelengths, we have observed dis- 
crepancies between the altimetric geolds and gravimetric geolds. While 
we could not Isolate the cause of this difference, It Is likely to be 
partially due to radial orbital uncertainty as was demonstrated for 
SKYLAB by Marsh et al. (1976). Any such errors restrict the utility of 
satellite altimetry In obtaining Information about the longest wave- 
lengths of the earth's gravitational field. At the present time ground 
tracking of orbiting satellites has provided the best Information on 
these longest wavelengths. In general, however, this is restricted to 
wavelengths with degree and order less than approximately 10. According 
to Khan (1976) the estimated accuracy In determination of Individual 
spherical harmonic coefficients for degree 11 and higher Is less than 
60X. More recent gravitational field solutions such as 6EM-9 (Lerch et 
al, 1977) ifiay however be more accurate. 

In order to answer the following question a hypothetical experiment 
can be performed. Compare an altimeter determining the geoid with a 
gravimeter measuring gravity. How well does each instrument perform in 
detecting Its respective signal? The answer to this can be found by 
examining the relative properties of the geoid and gravity field and the 
comparative accuracies of altimeters and sea gravimeters. 

Consider hypothetical measurements of geoid and gravity on a plane 
two dimensional earth. On this idealized body the geoid height is 
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estimated by utilizing sea surface measurements from a radar altimeter. 
Due to the presence of oceanographic noise we assume uncertainties in 
measurements of about 1 meter. Gravity is measured by a sea gravimeter 
to an accuracy of 10 mgal limited by cross coupling and other errors. 
Those are both somewhat arbitrary and conservative estimates but hope- 
fully reflect the measurement situation to a certain degree. While 
radar altimeters can measure sea surface height tr nossibly 10 cm 
(SEASAT-A), due to unknown oceanographic noise the geoid can only be 
determined to about 1 meter accuracy at present. Likewise, with marine 
gravimeters a precision of several mgal is possible, but in analyzing 
data due to two different ships at crossover, errors of 10 mgal are 
common. In our hypothetical experiment these instruments detect the 
geoid and gravity on a 5 km deep ocean overlying a 6 km thick crust. 
There is topography compensated in an Airy manner; the amplitude spec- 
trum of the topography is white with an amplitude of 3 km. The noise 
spectra for each instrument is also considered to be white. 

Then on the surface of the water the amplitude specti of gravity 
is obtained from equation 2 and spectrum of geoid is obtained from 
equation 3. A plot of these functions is in Figure 10. The vertical 
scale has been adjusted so that 10 mgals is equivalent in height to 1 
meter. Thus whenever one function is plotted higher than the other, it 
indicates the higher function has a larger signal to noise ratuc. From 
this plot we see that for wavelengths longer than 628 km, an altimeter 
will have a higher signal to noise ratio. For wavelengths shorter than 
628 km a marine gravimeter will record with a better signal to noise 
ratio. What we wish to indicate by this analysis is that below a 
certain wavelength a marine gravimeter can determine the gravitational 
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Comparative accuracies of an altimeter and gravimeter In 
detecting the gravitational field on a plane earth. This 
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field with a better signal to noise ratio than a radar altimeter. For 
longer wavelengths a radar altimeter performs better. This critical 
wavelength depends on the noise levels of the individual instruments, in 
our hypothetical example it was 628 km. It should be understood though 
that a radar altimeter still records short wavelength information; this 
was demonstrated in our previous example over the Aleutian trench. As 
radar altimeters and marine gravimeters improve in their ability to 
detect the gravitational field, there will be a different critical 
wavelength below which a gravimeter performs better. 

SUMMARY AND CONCLUSIONS 

In comparing the GEOS-3 geoid estimates with the gravimetric 
geoids over the different oceans, several features emerge. Constant 
offsets and tilts cause large-scale discrepancies of up to 24 m rms 
(Table 1). These errors can be due to radial orbital uncertainty, scale 
errors in the geoid, lack of atmospheric correction in constructing the 
geoid, or bias errors in altitude determination by the altimeter. We 
did not discern the relative importance of these factors in causing the 
discrepancies. Over large scale topographic features such as the 
Aleutian Trench there can be a large geoid signal with wavelengths 
shorter than can be resolved by 1® X 1° gravimetric geoids. 

In regions where there is good coverage of marine gravity data, and 
sufficient energy in the geoid at shorter wavelengths (less than 200 km) 
more detailed gravimetric geoids should be utilized. Examination of the 
topography and techniques utilizing the admittance between geoid and 
bathymetry should be useful in estimating the geoid signal at these 
short wavelengths. 
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After elimination of long wavelength differences, adjusted GEOS-3 
data still show discrepancies with the gravimetric geoids. Over the 
trench systems of the western Pacific there were differences of up to 
10 m. This is due to averaging gravity data over too large an area. In 
other regions there are discrepancies of several meters. Such dif- 
ferences indicate the need for estimation of the errors in geoid com- 
putations and the best possible orbital determinations for the GEOS-3 
satellite. 

Although some difficulties exist with constant offset d' q 
wavelength discrepancies, the GEOS-3 radar altimeter appear- .^ect 
geological features such as deep sea trenches and is an excellent instru- 
ment for acquiring measurements of the shape of the ocean surface. With 
the rapid mapping of the marine geoid on a global basis by GEOS-3 there 
will be a wealth of new information about the earth's gravitational 
field. It should be understood though that in detecting very short 
wavelength features of the field, a marine gravimeter will perform 
better than a radar altimeter. This is due to the long wavelengths 
being predominant in the geoid, while short wavelength anomalies are 
most obvious in gravity measurements. 
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Chapter III 


Techniques for Interpretation of Geo id Anomalies 
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Introduction 


When a new scientific instrument is developed and utilized for the 
study of the earth, there are always new and exciting discoveries. Such 
an instrument is the GEOS-3 radar altimeter, which has been used to make 
measurements of the shape of the sea surface over most of the worlds 
oceans. After corrections for errors and oceanographic effects this 
data is a determination of the marine geo id. 

Given the shape of this gravitational equipotential surface inferences 
can be made about the density inhomogeneities within the earth which cause 
geoid undulations. Prior to satellite altimetry determination of the geoid 
was the primary goal. Now with geoid measurements via satellite altimetry 
new questions can be asked. Why does the geoid have its particular shape and 
what are the structures within the earth which cause geoid anomalies? In 
Older to utilize altimetry data in studying the earth's interior mathematical 
techniques are necessary. For this reason we have developed several new 
analytic techniques for interpretation of geoid anomalies. 

Our procedure is to initially compute formulas for the geoid anomaly over 
idealized bodies. These serve to demonstrate various properties of geoid 
anomalies. Then formulas are developed for computing the geoid anomaly over 
an arbitrary two dimensional body. Most general procedures require the calcu- 
lation of geoid anomalies over three dimensional bodies. In order to under- 
stand the relationship between the gravity anomaly and the geoid, techniques 
in the frequency domain are developed. Finally, as practical examples, 
calculations of the geoid anomaly across continental margins and over sea- 
mounts are compared with actual geoid measurements. 
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Basic Formulas and Ideal Bodies 


In dealing with geold anomalies we are concerned with the computation 
of the shape of an equipotentlal surface, this surface being primarily 
ellipsoidal with small undulations due to an anomalous potential. 

Anomalous potential Is the difference between total potential and the 
normal potential due to the reference ellipsoid. If we know the anomalous 
potential, then the geold anomaly Is given by Brun's formula (Helskanen 
& Moritz. 1967). 


N ■ T 

7 ( 1 ) 

Where N Is geold height, T Is the anomalous potential, and y Is normal 
2 

gravity -980 cm/sec . A point mass will have an anomalous potential of 
G M 

T ■nr (2) 

-8 3 2 ' ^ 

where G Is the Newtonian constant 6.67 10 cm /gm sec , M Is the mass, and 
R Is the distance between the point mass and the observation point. This 
expression has two conventions: potential is zero at Infinity and positive 
everywhere else. Thus an excess mass corresponds to positive potential and 
positive geold height. Our basic problem then Is how to compute the geold 
anomaly due to an assemblage of point masses. 

i r / y* p d X dy dz iz 

N(x',y',z') - Y /// f- ,2,2 'Sr 

<J'JJ |(x-x') + (y-y')^ + (z-z')J 

Primed coordinates Indicate the observation point, unprimed coordinates are 
the integration variable, and p Is density. Initially we canpute the geold 
over ideal bodies; lines, sheets, and rectangles in both two and three dimensions. 

A finite horizontal line located at y^, with linear density \ (qm/cm) 
has a corresponding geoid anomaly of (Figure la) 
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6 
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X dx 


‘1 


[(x-x')2 + (yg-y')^ + (zo-z')^]^ 


( 4 ) 


This, when integrated, becomes 
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( 5 ) 


ga 

N > — log^ 
Y 


^(*2"*') ♦ + (yo-y')^ 

(XpX') + |(X,-X')2 + {yQ.y')2 



Figure IB shows a plot of this function, the geold height Is always positive 
and decays to 0 In the limit as |y| goes to Infinity. 

If the finite horizontal line Is Integrated along the z axis, one obtains 
the geold anomaly over a finite vertical sheet (Figure 2a). This Integral Is 
(see appendix for details) 


Go 

N - — { (X 2 . 22 ) - fil (*2.2]) - fi] {*1.22) * S] (x ,2 ) } 

Y • ' ' 


( 6 ) 


where tf Is the surface density (gm/cn?) and 


Gi (x,z) ■ (2-2') log^f(x-x') + |(x-x')^ + (yo-y*)^ } 

u j7j 

(2-2') ♦ (x-x') logg^(2-2') + [(x-x')^ + (yo-y')^ } 

- (X“X')^ - (yQ-y')^ - (x-x') |}x-x')2 ♦ (yg-y')^ *( 2 - 2 ')^ 
j(x-x') + |(x-x')2 ♦ (yo-y')2 +(z- 2 ')^ £x-x')2 + (yo-y*)^ 


sin 


and restricting ( 2 - 2 ') >0 

A plot of this function Is In Figure 2B. While always positive, the geold 
height decays at a slower rate than the finite horizontal line. 

Our next Ideal body Is an Infinite horizontal line, while the actual 
Integration Is quite simple we utilize a slightly longer derivation In order 
to show the relationship between gravity, deflection of the vertical, and the 
geold. This Is also the basis for a later derivation In the frequency domain. 
In the two dimensional geold calculations a complication arises, involving 
the zero reference level. Calculation of the geold anomaly over the Infinite 
horizontal line will Illustrate this complexity. 
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Calculation of geoid anomaly over a finite horizontal line and 
its resultant geoid height. 
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Calculation of geold anomaly over a vertical sheet. Resultant geold 

heights for the three dimensional (finite) and two dimensional 

(Infinite) vertical sheet are plotted. 

Figure 2 
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Due to an infinite horizontal line located at y»0, z*Zq. the horizontal 
conponent of gravity on the z*0 plane is 


■ -2 6 X v v 
g. 


( 8 ) 


(y‘ ^0 ) 

Deflection of the vertical is defined to be the angle between a vertical 
line along the z-axis and the local direction of gravity (Figure 3). On a 
plane earth the gravitational attraction has a vertical component of 
Y + Fg'^, where y is normal gravity and Fg'' is the vertical component of the 
gravitational attraction due to anomalous mass. The horizontal component is 


simply Fg", thus 


tan 6 * 


Fg' 


( 9 ) 


Y + Fg'^ 

v^ere 6 is deflection of vertical. Because Fg'^ « y. and 6 » 0 

Fg’^ -2GX y 

6 • ■ — 

Y Y (y^ + 

Deflection of the vertical in the y direction is just the slope of the geoid. 

a N -2GX y 

Y (y^ + Zq^) 


( 10 ) 




( 11 ) 


3 y 

Solving for geoid height N we obtain: 

G X 


N 


lo9g (y^ + Zg^) 


y 

yo 


(12) 


where y^ is an arbitrary constant of integration. With the two dimensional 
potential the usual procedure is to let the potential be zero at unit distance, 
or when + Zq^) » 1; it must be zero at some finite distance. This has 

two effects; it causes the absolute amplitude of N to depend upon the units 
used, and it Introduces both positive and negative geoid iieights for a positive 
mass. In dealing with two dimensional bodies it is necessary to ignore this 
arbitrary constant level. Thus in two dimensional models only 
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TWO DIMENSIONAL LINE 


“ g+F« 


Calculation of geoid anomaly over an infinite horizontal line. In 
B the horizontal attraction of gravity in F the vertical attraction 

g 

due to the line is Fg'^. Acceleration due to the earth is g, and 
deflection of the vertical is d. 

Figure 3 
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deviations from a constant level are important. 

After setting the potential at zero for unit distance we obtain for the geoid: 


G A 

logg (y2 + z^) 

Y ^ 


(13) 


When the infinite horizontal line is integrated along the z axis, one obtains 
the geoid anomaly over a two dimensional vertical sheet. 

G o loge [(yo-y'}^ + (z-z')^] dz 


N = 

and after integration 


G o 
yJz•^ 


(14) 


N = - {(z-z‘) logg [(yn-y')^ + (z-z')2] 


- 2 z + 2 yo-y'l Tan’^ 


.1 iz-z') 


\yo'y‘ 


(15) 


This anomaly is shown in figure 2c, for comparison purposes in figure 2b there 
is a geoid anomaly over a finite vertical sheet. Obviously they are quite 
similar in shape. This shows that two dimensional methods might approximate 

three dimensional ones in spite of the arbitrary height difference. 

A two dimensional rectangle is obtained by an integration of the mass 
line over both the y and z axis (Figure 4a). 


N 


^ r^2 ryz ^ [(y-y')2 + (z-z')2] dy dz 


^1 J yi 

after integration (see appendix for details) 

G a 


(16) 


N = [G 2 (y 2 » 22) - G 2 (yi, Z 2 ) - G 2 (y 2 , zi) (17) 

Y 


+ G2 (yi , zi )] 


u 



z 


Calculation of geoid anomaly over a two dimensional rectangle. 
Note that the geoid anomaly is never flat, rather it is always 
curving. 

Figure 4 
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where 


( 18 ) 


G2 (y. 2) » (y-y') {(z>2') logg [(y-y‘)^ + (z-z')^] - 2 z + 2 | y-y‘| 

(z-z') (y^y') 

Tan-1 } -2 y z + [(z-z')2 + (y-y')^] Tan"^ -+ (z-z‘) (y-y') 


|y-y'l 


(z-z‘) 


assuming {z-z') > 0 

In Figure 4b there is a ploi of this anomaly. Over the rectangle the 
geoid is not flat, it is always concave downwards until a point of inflection 
above the edge of the body. The implication of this curvacure is that in 
constructing geoid models the anomalous mass must integrate to zero, otherwise 
there will be a long wavelength curvature of the geoid anomaly. To make the 
anomalous mass sum to zero, both negative and positive densities are utilized. 


Integral Formulas For Two Dimensional Bodies 

While analytic solutions for ideal bodies have an initial utility in under- 
standing properties of geoid anomalies and verifying other formulas they are 
of limited utility in learning about the earth. Better models of the earth 
involve specifying densities of bodies of arbitrary shape and computing the 
resultant geoid anomaly. In order to achieve this result we have developed 
several techniques for the computation of geoid anomalies over bodies which 
are polygons in either two or three dimensions. These enable the models of 
virtually any shape to be calculated, the desired accuracy in shape only 
limited by the number of specified sides. 

A simple numerical solution to computing geoid anomalies was originally 
developed by Taiwan! , Poppe and Rabinowitz (1972) in order to compute the geoid 
across a two dimensional structure from the gravity. For a given polygon in 
two dimensions the horizontal component of gravity is computed with the line 
integral method of Taiwan! et al. (1959), then after division by normal gravity 
the deflection of vertical is obtained. 
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09 ) 


Fg" 

6 ■ 

Y 

As in our previous derivation for the infinite line, this is the slope of the 
geoid 

dn 

fi * 

dy 


( 20 ) 


so finally 


N = 


/: 


<5 (y) dy 


(21) 


This final integration is done with the trapezoidal rule or else Simpson's 
formula. Although it is a numerical solution it produces very accurate results 
in excellent agreement with exact analytic solutions. 

For geophysical interpretation of two dimensional bodies of arbitrary 
shape a method has been developed to directly compute the geoid anomaly over 
a polygon in two dimensions. In this technique the body is specified by the 
location of each vertex and the density of the body, as in Figure 5a. For 
any arbitrary two dimensional body the geoid anomaly is 


N = 


■‘iff 


logp [(y-y')^ + {z-z'fy dy dz 


( 22 ) 


To obtain a direct analytic solution for a two dimensional polygon we utilize 
Stokes' theorem to convert this surface integral to the following line integral; 


N 


-if 


{(y-y') logo [(y-y')^ + (z-z')^] - 2 (y-y ) + 2 (z-z') 


Tan 


-1 


(y-y’) 


( 23 ) 


} d z 


(z-z') 

and restricting (z-z') > 0 
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EEPRODUl i;;! Y OF THE 
ORIGINAL 



POLYGON 


with a polygonal approximation each side Is defined by the formula: 
z ■ m^ y + b^ 

then the contribution to the geold anomaly due to the ^th line segment Is: 
(see appendix for details) 


hi ■ (y-y')^ loQe C(y-y')^ + (z-z')^] - m| (y-y')^ 

+ (y-y') - ">1 (mi^-f) A 2 ^ log^ C(y-y')^ + (z-z')^] 




2(1-Hn|^) 


2x2 


j~ 0 ^ 1 ^) (y-y') + m^A 2 j 


-2 5)1^ A 2 ^ Tan P (y-y') 

(1-Hn^2)2 


-m*(y-y')2 + i-j^-l (y.y ■ ) 


-C^ 0^2 ]og^ 

(1+c^2)2 


p (1 + Ci^) + 2 Ci Di Z + Di^ "I 

L r“-i 


+ + (l-Ci^) Di2 Tan’Vyiyil^ 

(1+Ci2) (l+c^2)2 V z / 


yi-i 


(24) 


with A2 ■ y' +• -z' 

c^ * 1 / 

Di » -Bi/mi -y' 

" slope of line segment 
Bi = Intercept of line segment 
Z = m^ y + bi 
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(25) 

L * number of segments 


thus the total geold anomaly 1s 



y 1-1 

Due to singularities In the Integral formula the following special cases need 


also be considered: 

a) m^ ■ 0 then n^ ■ 0 (slope of line segment Is zero) 

b) C{ ■ 0 then (slope of line segment Is Infinite) 

■ {(y^-y‘) C( 2 -z') loge C(yi-y')^ + (z-z')^3 - 2 (z-z') 

+ 2|y,-y| -Myi-y ')2 


|(yi-y')| 

+ [(yi-y*)^ + (z-z')2] Tan*^ 


(yi-y') 

(z-z') 


+ (yi-y') (z-z')} 


c) B^/m^- + y' -0 


zi 

zi-1 


(26) 


then 


(the extension of line segment passes 
through observation point) 


3z2 


"1 “ ^°9e C(l+-— t) z 2] + z^ Tan"^ ( )} 


1 


2mi 


mi 


2mi mi 

assuming z' ■ 0, z > 0 


(27) 


For a given geophysical model a number of Individual bodies are specified 
by their respective densities and shape as defined by a number of vertices 
of a polygon. For each body the geold effect Is calculated using either the 
numerical technique In formula (21) or the analytic technique of formula (24). 
After summation over all bodies the total geold height due to the model Is 
obtained. In practice both of these mathematical techniques yield Identical 
results. In Figure 5b there Is a sample computation of geold height and 
vertical component of gravitational attraction over a pentagonal body. The 
geold has a shape similar to the gravity anomaly but Is much wider; the 
geold height decays with distance at a ^ower rate. 


In modeling two precautions are necessary. First the constant level of the 
geoid height must be ignored. This is because of the arbitrary location of 
zero potential for a two dimensional body. A second precaution concerns the 
specification of densities. In geoid calculations the anomalous density must 
be utilized rather than the total density. Hiis is because Brun's formula 
(formula (1)) requires anomalous potential rather than total potential. If 
a given model is in isostatic equilibrium this second precaution requires 
that the sum of the anomalous masses In a given column be zero. This require- 
ment is different than in standard gravity models If the total mass did not 
sum to zero then a constant gravity value could always be subtracted from calculated 
gravity to obtain the gravity anomaly. In geoid computations this procedure 
does not work, the reason being that a two dimensional rectangle (as discussed 
previously) does not have a flat geoid anomaly above it, whereas in gravity 
calculations a two dimensional rectangle, except for edge effects, has a flat 
gravity anomaly. 

As a practical example of the utility of these methods we consider the 
study of an Atlantic type continental margin. In the transition between the 
deep sea and the continental mar jin of Nova Scotia In Eastern Canada, the 
ocean bottom changes in depth by five kilometers. Likewise the 
ocean surface changes in height by six meters (Figure 6). This change in 
geoid height has been determined from the GEOS-3 satellite. A natural question 
arises then, what is the cause of this change in geoid height? To help answer 
this question, a two dimensional density model of the Nova Scotia margin has 
been constructed (Figure 6). This is a simple three component model with water, 
crust, and mantle material. It is fully compensated in an Airy manner. Util- 
izing the two dimensional technique described in this section, we have calcu- 
lated the geoid anomaly over this model. As Figure 6 indicates, the calculated 
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Figure 6 

Geo1d anomaly across Nova Scotia margin. 
Is at 44.5°N, 63.25°W and right end Is at 
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500 km 


WATER ./»=I.O 



t end of GEOS- 3 profile 
i“N, 58.7“W. 




geold height Is In good agreement with the observed GEOS-3 data. It can be 
concluded that a simple density model of the Nova Scotia margin which Is In 
Isostatic equilibrium Is capable of successfully modeling the observed GEOS-3 
data. Due to the noise level of the data a more elaborate model was believed 
unjustified. 
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Geold Calculations For Three Ditnensional Bodies 


While the two dimensional approximation has a certain utility, more 
realistic geophysical modeling requires th<t computation of the geold anomaly 
due to bodies of arbitrary siiape and density In three dimensions. In the 
method now developed the shape of the body Is prescribed to be a polyhedron. 
Any shape can be approximated by this method by simply Including sufficient 
detail with many facets. This same technique was developed in gravity 
modeling by Coggon (1976), Paul (1974), and Barnett (1976). 

In Figure 7a we Illustrate a polyhedron which Is assumed to have constant 
density. To calculate the geold height for this body the fundamental Integral 
formula 3 has to be solved. Conceptually our method of solution Is quite 
simple. This volwne Integral Is converted to a surface Integral by the 
divergence theorem. For each facet of the polyhedron this surface Integral 
can be solved by conversion to a line Integral via Stokes' theorem. Thus 
the total volume Integral Is calculated by a summation of surface integrals, 
one for each polygonal facet. Each surface Integral is calculated by a 
summation of line integrals, one for each line segment on a given facet. 

This Is an exact analytic technique and Is accurate for any shape which can 
be approximated as a polyhedron. 

Rewriting equation 3 In the form 

where r is the unit vector In the 
radial direction 
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THREE DIMENSIONAL POLYHEDRON 



Calculation of 9«o1d anonaly over a thr« dimensional polyhedron. 

Figure 7 
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(>!■' 


TW-: 


I 


( X - x') X 


>') V * [z-i'AA 


C(x - x ')2 + {y - y')^ + (2 - 2')^]* 
and 7 1 s the vector del operator 

We see that It can be changed to a surface integral with the divergence theorem. 


Thus: 


N ■ G oj^ r • U 4 ds 


O'" " G p jh^ C(x - X') + (y - y') U42 ♦ (2-2') U43] ds 

s [(X - x *)2 + (y - y')^ + (2 - 2 ‘)¥ 

where u^ Is a unit vector normal to the element of surface area ds, and has 
components (u^^ , u^g* “ 43 )- 

This equation now Is a surface Integral and mu*, je solved for each 
Individual facet upon which the unit normal vector U 4 Is constant. After this 
Integral Is computed for each face, the results are summed for all facets 
thereby calculating the geold anomaly due to the entire body. Thus 


N » r I 
1 1 


where Is equation 29 evaluated for 
the 1th facet of the polyhedron. (30) 


Solution of this Integral for each face of the polyhedron requires a coordinate 
transformation, rotating the facet until It is normal to cne of the 
new axes. This technique was originally used for solving integrals in gravity 
and magnetic computations by Barnett (1976) and also Paul (’974). 

This coordinate transformation Involves a rotation and translation. To 
construct the necessary rotation matrix, three unit vectors must be first cal- 
culated. For a giver facet (Figure 7b) these three vectors are determined In 
the following way. Each vertex In the facet is numbered in counter clockwise 
order when looking Into the body. Unit vector u^ extends from point 1 to p^int 
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2 , and unit vector U2 extends from point 2 to point 3 . Given this numbering 
scheme u^ x U2 (vector cross product) Is the outward facing normal for the 
facet. With these vectors, two additional unit vectors can be calculated 


■ u^ X U2 


U3 - u^ X ui 


11 

U12 

“13 

'31 

U32 

“33 

41 

“42 

“43 


where u^ Is the 
outward unit normal 

An orthogonal rotation matrix Is constructed from these unit vectors. 

|u, 
u 
u 

n I ■»*- T, 

where u^j Is the jth component 
of the unit vector u.{. 

A coordinate translation Is also applied so that the observation point (x', y' , 2' 
becomes located at the new zero origin. The complete coordinate transformation 
changes the (x, y, z) coordinates to a new (c, n> a) system. In which 
(x', y', z') becomes (0, 0, 0,). 

Thus 


(31) 


where 

c' » u^i x' + u^2 y' + ui3 z' 

n' • U31 x' + U32 y' + U33 z' 

0' * u^T x' + U42 y' + U43 z' 

After coordinate transformation the Integral 29 for the 1 'th facet changes 
to the simpler form: 


c 


“11 

“12 

“13 


X 



n 


“31 

“32 

“33 


y 

- 

n' 

e 


%1 

“42 

“43 


z 


0' 


ii ’Ll 


// 


1 ef d ; d »i 

I 


+ 0^^ 


(32) 
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In order to solve this surface Integral for the 1'th face, we change it to 

a line Integral over the sides of the polygonal facet via Stokes' theorem. 

•4 


I 




( 33 ) 


where the Integration Is counter 
clockwise (when looking Into the body) 
around the line segments bounding the 
rth face of the polyhedron. It must be 
a closed contour. 

Each side of the facet Is a line segment with an equation 

Z * mjfi + (34) 

where m^ Is the slope of the j'th line segment 

"*j “ ^^j+i ■ 


'’j+1 


- H: 


and also = c . - m. n. 

J J J J 

Thus the total Integral for the 1'th facet Is a sum of line integrals over 

all of the bounding line segments. 

1 

j 




(35) 


and 


^j+1 

* G_e. j ®i ^°9e 


(36) 




(jmjn + Cj)^ + ^dn 

and Lj Is the line Integral for the j'th line segment of the 1'th facet of 
the polyhedron. 

This assumes that m^ Is finite . If It Is not then dn is zero and consequently 
Si s Lj. Also if 0 -j is zero, then I'j is zero. After integration of equation 
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36 and summation over all line segments of all facets, the final result Is 
obtained. (See appendix). 

N - 01 ^ 

1 ^ j 


loQg ^(m^r 




Omjn + + 



fs™, 







( 37 ) 


+ n 





(n!j0^ - Cj»i) 
jjmjT) + Cj)^ + 



’’j+1 


n 


j 


This expression then. In conjunction with the coordinate transformation 
of equation 31, when evaluated, gives the geold anomaly over an arbitrary 
polyhedron. As a practical matter. It Is difficult to specify the coordin- 
ates of the body In the orderly manner required. To simplify matters all 
faces could be triangles and the body would be assembled In the manner des- 
cribed by Barnett (1976). Also It should be noted that both the two dimen- 
sional and three dimensional solutions to geold anomaly modeling are exact 
analytic expressions. This Is In contrast to the asymptotic expansion 
technique for geold anomaly modeling developed by Ockendon and Turcotte (1977), 
which although easier to evaluate requires the body to be very thin. Another 
modeling technique was utilized In a study of a seamount region by Fisher (1976). 
In many studies of geold anomalies It will be necessary to compute the resultant 
anomaly on a spherical earth. The three dimensional method developed here 
can be utilized for that purpose, as long as it Is realized that the connecting 
segments between points on the polygon are linear. An alternative approach 
would be to utilize the method developed by Johnson and Litehlser (1972), 

who give an expression for the potential due to a body on a spherical earth. 
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Their method In combination with Brun's formula would be a good numerical 
solution to geold modeling on a spherical earth. Another numerical technique 
for computing geold anomalies on a spherical earth was developed and applied 
In a study by Bowin et. a1. (1975). 

As a geophysical example, this three dimenclonal modeling technique Is 
utilized to study a geold anomaly which has been observed above the Truk Islands 
(part of the Caroline Islands In the western Pacific). This geold anomaly Is 
5 meters In amplitude and Is above a seamount which Is 4 km high . Quite 
simply this geold high can be explained as due to the excess mass of t1« 
seamount Itself, and the deficit mass of Its compensating body. To model this 
anomaly the topography was digitized and a polyhedral body was constructed* 

In addition a polyhedral compensating mass was calculated. An Airy type crust 
with a compensation depth of 30 km was utilized. In Figure 8 a cross section 
of this model Is depicted, along with both the calculated and observed geold 
height. From the good agreement of the two. It can be said that this geold 
anomaly Is explained by an Isostatic model of this seamount. 
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sediment and 6 km basement. 


Techniques In the Frequency Dwialn 

A different understanding of the geold Is possible In the wavenumber 
domain. This Is done by examining the relationship between the Fourier 
transform of gravity and the Fourier transform of the geold. In other words, 
what Is the transfer function between gravity and the geold? Our derivation 
Is analagous to the derivation of formula 21, but Is done In the frequency 
domain. 

Given the vertical component of gravity due to a two dimensional body 
we go to the wavenumber domain via the Fourier transfonn. 

+«• 

flf/) ■ f f/ (y) -y (38) 


To obtain the Fourier transform of the horizontal component of gravity we 


Hilbert transform 


9 

■ 0 ^) 




(39) 


where multiplication by 1 sgn (k) Is 
the Hilbert transformation (Papoulls, 1962) 
also sgn (k) = +1 k> o and 1 « £- 13 ^ 

-1 k< 0 


This procedure can be verified by comparing the Fourier transforms of the 
vertical and horizontal components of gravity due to an Infinite mass line. 

After division by normal gravity we get the Fourier transform of the deflection 
of the vertical. 

. ) sgn (k) 

Y ~ 

Utilizing the Integration theorem In the frequency domain the Fourier transform 
of the geold Is obtained. This Is because the geold Is the Integral of deflection 
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of the vertical . 


^(N) • ») . 1 san (k) 

3 1 Y k (41) 

or ^CN(y)] • 1 ^ [F ’ (y)] 

yTm ^ 

Thus given the Fourier transform of gravity we multiply by a transfer function 
S(k ) to obtain the Fourier transform of the geold. For the two dimensional 
case: 


S(k ) 


TTH 


(42) 


This transfer function shows that a geold anomaly over a body will have 
the same phase spectra as the gravity anomaly. In simpler words this is our 
previous observation. A geold anomaly over a body has the same general shape. 
Any maxImiBii In the geold occurs exactly where a gravity maximum does. However 
the geold anomaly for a given wavelength Is changed In amplitude by the factor 
( 1/1 ItD- Thus for long wavelengths ( k«l ) the geold anomaly^will be 
magnified, while at short wavelengths ( k»l ) diminished. This Is our 
observation that geold anomalies are "wider". 
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Conclusions 


In this paper we have been concerned with developing the techniques 
necessary for the Interpretation of geold anomalies. These are mathematical 
formulas which for a given geological model can be utilized to compute the 
resultant geold anomaly. If the model Is considered two dimensional then 
densities for a polygon In two dimensions are specified and using formulas 
24 and 25 the geold heights are calculated. In this method only the relative 
heights are Important. Any geologic body can be modeled with a polyhedral 
shape In three dimensions. If density Is variable within the body several 
polyhedra may be utilized. After specifying densities for this, formula 37 
Is utilized to calculate the resultant geold anomaly. To understand geold 
anomalies, It Is seen that they basically have, the same shape as gravity 
anomalies, but with the long wavelength components amplified. Techniques 
"such as these are only tools. Their real Importance Is in enabling the 
Interpretation of observed geold anomalies, such as are detected by GEOS-3. 
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Geold Anomalies Over Oceanic Features 
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Geold Anomalies over Oceanic Features 


In using satellite altimeter measurements to study the Earth's Interior 
a simple question can be asked. What Is the shape of the ocean surface over 
bathymetric features? Physically such features would be expected to have 
related geold anomalies because there Is a large density contrast between 
oceanic crust and the surrounding water. At the bottom of the ocean floor 
there are certain types of large structures; seamounts, continental margins, 
fracture zones, mid-ocean ridges, and deep sea trenches. Over each one of 
these features there are characteristic undulations of the ocean surface. 
Being directly related to the density structure, these measured shapes can 
be utilized to infer the internal lateral density variations within the 
Earth. Using the GEOS-3 radar altimeter geold anomalies have been identified 
* over these structures. 

Bordering the continents in the Atlantic ocean are passive continental 
margins. Detailed geophysical studies of these regions usually involve the 
study of several Important questions; location of the ocean-continent 
boundary, nature of the mechanism of isostasy, and history of subsidence 
along the margin. Of these questions, geold measurements may be of 
assistance with the first two. As an example of a passive margin consider 
the East Coast United States. Across this margin the geold has a 6 meter 
step similar In sense to bottom topography (Figure 1). In plac’S this 
anomaly Is quite distinct, as off Nova Scotia; further south the increase 
is quite gradual. An important area of research would be to ascertain why 
this shape varies along the coast. This study would probably be concerned 
with isostatic geoid anomalies. 
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At many of the coi *1nental margins In the Pacific Ocean active 
subductlon of the lithosphere occurs at deep sea trenches. These are 
active continental margins. At these margins there Is a characteristic 
geold anomaly. With GEOS-3 data we see that above the trench axis the 
ocean surface forms deep valleys up to -20 m In amplitude. Over the Aleutian 
Trench the sea surface has a depression of -13 m, while at the Mariana 
trench the geold dips by almost -20 m (Figure 2). Seaward of the Aleutian 
Trench there Is a long wavelength Increase In the geold towards the Aleutians. 
This anomaly Is concave downward and continues to Increase landward of the 
trench where the depths are shallower. Because It Is the nature of the 
geold to show long wavelength anomalies It Is expected to be most sensitive 
to deep density Inhomogeneities rather than crustal structures. As the 
downgoing slab Is a large and deep Inhomogeneity, any anomaly caused by 
It would be more evident In the geold rather than in gravity. 

Another primary type of geold anomaly is above mid-ocean ridges 
where accretion of the lithosphere occurs. One example of this type 
anomaly Is across the Southwest Indian Ridge in Figure 3. Gravity and 
topography profiles are from a nearby ship track; the geold data are from 
the GEOS-3 satellite. Although gravity data show little of a long wavelength 
anomaly, the geold has a distinct anomaly centered over the ridge axis. 

This Is about 7 meters In amplitude and at least 1000 km wide. A track 
chart of altimeter data acquired over both the Southwest Indian and Southeast 
Indian Ridge is illustrated in Figure 4. Other geoid profiles across the 
Southwest Indian Ridge are shown in Figure 5. Each of these six profiles 
shows a higher region to the sout‘.east of the ridge axis followed by a sudden 
change in slope directly at the axis. As there appears to be a broad 
regional slope to the data a linear trend has been removed from each profile. 
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Geoid profiles over trenches and a marginal sea. 
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Profiles over the Southwest Indian Ridge. 
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Figure 5 

Geoid profiles across the Southwest Indian Ridge. Bathymetry 

and gravity profiles from Conrad cruise 11-04 are also shown. 
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These profiles are 1n Figure 6. Each profile over the Southwest Indian 
Ridge after removal of a linear trend shows a characteristic geold anomaly 
similar to the one In Figure 4. This anomaly Is maximum at the ridge axis 
and almost 10 m In amplitude. It Is 1800 km wide. Profiles across the 
Southeast Indian Ridge show a much larger regional background as shown In 
Figure 7. After removal of a linear trend the characteristic ridge anomaly 
Is evident In Figure 8. This anomaly has a 6 meter amplitude and 2000 km 
width. Presumably this ridge anomaly Is the result of density anomalies 
due to creation and evolution of the lithosphere. A detailed study of 
this and other ridge anomalies should be undertaken to ascertain any 
dependence upon spreading rate and any possible asymmetry. Using a thermal 
model of the lithosphere, Haxby and Turcotte (1978) have modeled the ridge 
anomaly In the North Atlantic. 

Other smaller oceanic structures have shorter wavelength geold 
anomalies. Above Islands there Is a geold anomaly several meters high. 

For example 6 meters over Bermuda (Figure 1) and 4 meters over the Caroline 
Islands (Figure 9). At fracture zones where the crust has a significant 
age difference on either side there will be a difference In depth which 
causes a geold anomaly. Over the Romanche Fracture Zone the geold has 
a step of 4 or 5 meters (Figure 10). With similar large fracture zones It 
will be possible to use geold height measurements to trace fracture zones 
where ship measurements are scarce. Another type geold anomaly occurs over 
the Bismarck Sea. This Is a small marginal basin several kilometers shallower 
than adjacent ocean. Above this region the ocean surface forms a broad high, 
450 km wide and 7 m high (Figure 2). 
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Geold profiles across the Southwest Indian Ridge with a linear 
trend removed. 
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Geoid profiles across the Southeast Indian Ridge. Bathymetry 
and gravity profiles from Conrad cruise 11-05 are also shown. 
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trend removed. 
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CAROLINE 

Caroline Islands geold anomaly. This anomaly over the Islands 

ISLANDS Is 3 or 4 m In amplitude, also there Is a step In geold height 


of 7 m across the Islands. 
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Chapter V 


Geold Corrections In Marine Gravimetry 
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Introduction 


Due to undulations of the geoid free-air gravity anomalies measured 
at sea are caused both by lateral changes in mass distribution and also 
by changes in height above or below the reference ellipsoid. In those 
gravity modeling studies which do not consider the latter effect of height 
changes, there may sometimes be serious misinterpretation. The required 
correction term to eliminate this effect is at most ±27 mgal and most 
significant when studying the geophysical implications of the long wave- 
length gravitational field. 

To understand the need for this correction, consider the standard 
procedure in gravity model inc; a density model of the anomalous mass is 
used to calculate the vertical component of gravitational attraction 
which is then compared with the free-air gravity anomaly. This procedure 
implicitly assumes the equivalence between free-air gravity anomalies and the 
gravitational attraction (in the vertical direction) of the anomalous mass. 
This assumption is only an approximation; due to geoid undulations there 
can be substantial errors in this procedure. This approximation is well 
known to geodesists (Bomford, 1962) but is not commonly acknowledged in 
geophysical studies. 

For models of localized mass anomalies or short wavelength features 
this approximation is probably valid. For regional studies, however, with 
the objective of interpreting longer wavelength components of the gravity 
field the importance of correcting the free-air gravity anomaly for geoidal 
undulations is greatly increased. 

For example, in a study of the North Atlantic regional gravity field 
the effect of this approximation was not discussed by Cochran and Talwani 
(1978). In this work there would have been a correction term of up to 
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approximately 20 mgal added to their observed free-air gravity prior to 
Interpretation. Other studies In the Japan sea and Izu-Bonin region 
(Segawa and Tomoda, 1976) and the Aleutian arc (Grow, 1973) have similarly 
not considered this additional correction. A modeling study of mid-ocean 
ridges (Lambeck, 1972) was one of very few to Include this Indirect effect 
In the theoretical calculations. Due to the extensive neglect of this 
effect In geophysical studies we briefly discuss Its origin and then Indicate 
the types of physical situations In which It Is most significant. 

Origin of Error 

When measurements of gravity are made aboard ship at sea they are of 
course actually being made on the surface of the geold. In some cases. 

In order to properly utilize this data for geophysical Interpretation, It 
Is necessary to correct for the undulations of the geold surface. Consider the 
standard method for obtaining the gravity anomaly. After measurement of 
gravity on the geoid, the normal gravity on the ellipsoid Is subtracted thus 
yielding the g**av1ty anomaly. 

A 9 “ 9obs ■ '^n where Ag Is free-air gravity anomaly 

gobs is measured gravity on the geold 
Yo Is gravitational attraction of the 

ellipsoid, on the surface of that ellipsoid 
It Is Important to note that for determining the free-air gravity anomaly, 
these two gravitational attractions are measured with respect to different 
references . The geold Is the reference for the observed gravity and the 
ellipsoid the reference for theoretical gravity. The difference In height 
between the ellipsoid and geold (Figure 1) contributes to the difference In 
gravitational attraction between what Is measured (g^jj^) and what is calculated 
(y^), and thus the resulting gravity anomaly Is not solely due to an anomalous 
mass. To see the need for this additional correction consider the following 
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Figure 1 


Measurements of gravity (9 qJjj) are made on the geoid which is at a 

height N above the reference ellipsoid. Attraction of the ellipsoid 

(y ) is calculated at a different height, on the surface of that 
n 

ellipsoid. This height difference is the reason for the geoid 
correction term, 
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expression which simply equates forces of gravitational attraction: 

ag • g . - y =^22. n ♦ gravitational effect of mass between geold 

obs n dr 

and ellipsoid + gravitational effect of anomalous mass 


I 

I 

V 

I 


where Is radial gradient of normal gravity 
N Is height of the geold above the 
reference ellipsoid - In meters. 

On the right side of this equation the first term Is a free-air correction, 
and the second term a Bouguer correction. The free-air term Is -.3086 mgal/ 
meter and using a planar approximation the Bouguer correction Is 
+.043 mgal/meter assuming a water density of 1.03 gm/cm3. Correction for 
the Bouguer term subtracts the effect of excess water mass for positive 
geold height and adds In a gravitational effect of the missing water mass 
for negative geold height. Thus to obtain the gravitational effect of 
the anomalous mass: 


gravitational effect of ancmialous mass (In mgal) = Ag + .2656*N 

This equation is obtained by transferring the free-air and Bouguer terms 
to the other side of the previous equation. The term which Is added to the 
gravity anomaly can be designated the geold correction term . This term 
1s numerically equal to .2656*N. This geold correction term is added to 
the free-air anomaly (Ag) In order to obtain the gravitational attraction 
of anomalous mass - In mgal. 

Significance of Correction 

This correction will be at most ±27 mgal, corresponding to the fact 
that geold height (N) Is rarely above or below 100 meters. To consider the 
importance of this correction it 1s necessary to understand the properties 
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of the geoid. Undulations of this surface are very smooth, with very long 
wavelengths. For many regions the corresponding geoid correction will only 
appear as a constant offset in the data, which in gravity studies is not 
always important. However, when considering the long wavelength gravitational 
field (wavelengths longer than approximately 500 kra) this correction becomes 
of greater significance in geophysical studies. At the present time, scientists 
are beginning to investigate the long wavelength components of the gravi- 
tational field, and it is in these studies that the geoid correction in 
marine gravity must be considered. 

To give an example of how this correction might be applied, we consider 
a profile of gravity across the Tonga- Kermadec trench (Figure 2). Above the 
bathymetry profile, the gravity anomaly (Ag) referred to the geoid is shown. 
Utilizing es. ’’’»tes of geoid height from the GEOS-3 satellite we have also 
illustrated the ^ old height along track. 

As the qeoic ,jrrection term is .2656*N it is directly proportional 
to geoid height (N). Thus the scale on the top left of Figure 2 is 

for geoid height (N) while the scale on the top right is for the geoid 

correction term which is a value in mgal to be added to the free-air 
anomaly (Ag). After addition of this correction to the gravity anomaly, the 
gravity referred to the ellipsoid is obtained and illustrated in the same 
figure. In this case the correction is at most 18 mgals. For this region we 

must affirm the existence of a regional gravity high behind the island arc. 

This is an important conclusion and directly shows the application of the 
geoid correction term . 

For the purposes of relatively short wavelength gravity modeling this 
correction may not be very significant while, for more regional studies, 
failure to take the geoid effect into account may lead to considerable 
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misinterpretation. The Importance of considering the geold height correction 
for measured gravity Is further emphasized as a better definition of the 
structure of the earth's crustal layers leads to efforts to determine the 
nature of deeper structures. Lateral density heterogeneities In the mantle, 
for example, produce gravity variations similar In wavelength to those of 
the undulations of the geold. As an example of the effect of the geold 
correction term on regional gravity and the consequences for Interpretation, 
one can examine the regional high associated with marginal basins behind 
western Pacific Island arcs. The regional free-air gravity anomaly field 
landward of the Mariana and Izu Bonin arcs averages approximately *ZS mgal 
while the Pacific ocean basin has a regional field of approximately zero. 

The geold height above the reference ellipsoid varies from as high as 60 meters 
behind these arcs to zero In the Pacific basin thus producing a geold 
correction term of up to 16 mgal. After addition of this term to the 
free-air anomaly the gravitational attraction of anomalous mass In the 
marginal basin Is 40 mgal, over the Pacific ocean basin It Is still zero. 

This would result in a different Interpretation about the extent and density 
of deep lateral mass heterogeneities. 

Summary 

The need for this geold correction term was stated previously by 
Bomford (1962, p 443), and Is thus not a new Idea. This Is only a timely 
reminder, since with the determinations of the geold by the GEOS-3 radar 

f altimeter the correction Is now feasible. Of course It was possible to 

i 

Implement this correction prior to radar altimetry by using global gravi- 
tational field models, however, with GEOS-3 geold data it can be done with 
! a higher degree of accuracy. In geophysical studies of the long wavelength 
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I gravitational field this correction Is mandatory, although due to the long 

I wavelength nature of the geold It may not be significant In studies of 

I the shorter wavelength gravitational field. 
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Introduction 

Measurement of the gravitational field In the world's oceans has 
resulted In an Increased understanding of the physical processes which 
occur In the vicinity of deep sea trenches. It has long been recognized 
that Island arcs and adjacent trenches are the site of some of the 
largest free>a1r gravity anomalies over the surface of the Earth. Ever 
since the development of pendulun measurements of gravity In submarines, 
a basic understanding cf the primary features of gravity at Island arcs 
has been possible. With only a few data points In a given transect It 
was possible to recognize the existence of a significant low In the free- 
air gravity anomaly directly above the trench axis and an adjacent high 
over the Island arc. Utilizing the Vening Melnesz pendulum apparatus on 
the submarine H.M.S. Telemachus . four crossings of the Tonga Trench 
Indicated this characteristic pattern of gravity anomalies (Talwani, et 
al., 1961). On one of the profiles a value of -224 mgal was observed 
over the Tonga Trench and a positive free-air anomaly of -•‘153 mgal over 
the Island arc. On the basis of earlier studies In the Indonesian archi- 
pelago Vening Melnesz (1954) observed that significant Isostatic gravity 
anomalies existed In the region and suggested that the trench system 
existed as a result of horizontally compressive and shear stress. 

When It became possible to make continuous gravity measurenents on 
board a moving surface ship, a wealth of new Information about the 
Earth's gravitational field became available. With this new data It 
became evident that many trenches had a distinct gravity anomaly seaward 
of the trench axis. This has been termed the Outer Gravity High (Watts 
and Talwani. 1974). One plausible explanation of this anomaly Is that 
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at a trench the lithosphere undergoes elastic deformation produc'M.? a 
bathymetric risw and consequent gravity high. 

With plate tectonic theory It became evident that as the litho- 
sphere descended Into the relatively hot asthenosphere. the slow rate of 
heat conduction permitted differences In temperature to be maintained. 

These thermal conditions result In a density contrast betvwen the cooler 
slab and the surrounding hotter mantle and a consequent gravity anomaly. 

It Is uncertain what the actual differences In density are, but the 
configuration of the descending slab Is known from seismic .cudles. 

Thus models of the density structure Including the downgoing lithosphere 
have been possible end were actually Cwr.itructed for the Chile Trench 
(Grow and Bowin, 1975), Aleutians (Grow, 1973), and Japan Trench (Segawa 
and Toffloda, 1976). In all of these studies each new development In 
Instrumentation or theory resulted In a better determination and explanation 
of the gravity field In the vicinity of deep sea trenches. 

Measurement of the gravitational field In a marine environment has 
heretofore been concerned with the radial component of the gravity 
vector. However, as the gravitational field is a conservative force 
field. It can equally be described In terms of the scalar gravitational 
potential. On the Earth the equipctentlal surface of the gravity field 
which coincides with mean sea level Is te>'med the geold. In the absence 
of noise, determination of the shape of the geold over the Earth's 
surface would give Information on the gravitational field equivalent to 
that obtained by gravity vector measurements. In the presence o^ noise 
geold height measurements give more useful Information for wavelengths 
between several hundred and several thousand kilometers (Chapman and 
Taiwan! , 1979). Determination of the shape of the geoid in marine areas 
can contribute new and valuable information about the gravitational 
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field and the Interior densities within the Earth. In this study 
measurvaents of geold height In the vicinity of deep sea trenches 
enable Identification of a characteristic geold anomaly. For the 
Tonga- Kermadec Trench, geold and gravity values are used to determine 
a new model of the density structure. These new measurements give 
better constraints on variations In mantle density than were previously 
available with only gravity measurements. 
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Data Acquisition 

Geold measurements are made by using a radar altimeter mounted 
on a spacecraft which orbits the Earth. In utilizing the GEOS-3 satellite 
to obtain estimates of geold height, a certain procedure Is used to 
collect and process the data. Mounted on board this orbiting spacecraft 
Is a radar altimeter which Is always oriented toward the center of mass 
of the Earth. During data acquisition the altimeter emits radar pulses 
which are reflected from the ocean surface and received by spacecraft 
antennas. In ch'^s ^^lanner measurements of the satellite altitude are 
I obtained 100 times per second. When combined with accurate orbital 

tracking of the satellite, the height of the Instantaneous sea surface 
relative to the reference ellipsoid Is determined. Due to the existence 
of forces other than from the Earth's gravitational field, the sea 
surface at any given time may depart from the actual geold. Such 
differences are caused by tides, currents, and other oceanographic 
phenomena. In general though, for the deep ocean these discrepancies 
can be considered noise which In amplitude Is less than approximately 
1 meter. After acquisition of the data, it Is time averaged over .1024 
second Intervals; a second smoothing procedure averages the data over 
2.2 or 3.3 second Intervals. The data which we utilize is an average 
value of geold height over a 14.4 km or 21.6 km distance; the ground 
speed of the sub-satellite track Is approximately 6.55 km/sec. For all 
of the data which we use In this study, a listing of the instrimient 
type, transmission mode, and date of acquisition Is available In Table 1. 
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TABLE I • DATA CHARACTERISTICS 
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t Orbital Accuracy 

A > lOi^ rw based on I day arc by Hal lops Flight Center - GEH-8 gravity nodel. 

0 • 3-IOu ras based on 1-8 revolution arc by Wallops Flight Center - GEM-B gravity 
J - 3u i«s based on 3 revolution arc by Naval Surface Weapons Center - NWL gravity 


GEOID ANOMALIES OVER DEEP SEA TRENC HES 

With the use of this GEOS-3 altimeter data, an extensive study of 
geold heights over Island arc regions has been undertaken. As the data 
was examined. It became possible to Identify a characteristic geold anom 
aly which exists over all major deep-sea trenches In the world's oceans. 
This geold anomaly Is caused by the lateral density differences within 
the Earth beneath the Island arc region. With aoequate additional 
constraints and proper Interpretation, It can provide useful Information 
about the density and structure of the upper mantle and the physical 
processes which occur at deep sea trenches. 

To discern the characteristic geold anomaly over deep sea trenches, 
we examine the various trench systems In the Pacific and Indian Oceans. 
These are In turn the Aleutian and Kuril Island arc system, the Mariana 
and Bonin Island arc-trench region, the Philippine Trench, the Tonga- 
Kermadec region, Peru-Chlle Trench, and lastly the Java Trench. Over 
each of these deep-sea trenches a set of representative geold profiles 
Is displayed. Track locations of all the data are illustrated in Figure 
1, while Instrument and data characteristics are listed in Table 1. For 
each of the geold tracks the data has been projected to a line per- 
pendicular to the strike of the trench ~.xls; as the original tracks are 
fairly linear, this only represents a scale change. Profiles from each 
trench region are displayed In Figures 2 through 7. Bathymetry and 
gravity data from a nearby ship track Is shown along with geold profiles 
from the GEOS-3 satellite. 

Geold profiles from the Aleutian and Kuril Trench are illustrated 
In Figure 2. These data clearly indicate the type of profile which will 
be seen as characteristic of all deeo sea trenches. Beginning about 
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2000 kilometers seaward of the trench axis, the geold has a constant 
Increase In height as the axis Is approached. This slope Is about 
20 meters over a distance of 1800 kilometers. Directly above the deep- 
sea trench there Is an approximate 13 meter decrease to form a low 
150 kilometers In width. Immediately landward of the minimum there Is 
an approximate 18 meter Increase until a maximum Is reached directly 
above the Island arc. Behind the arc there Is a gradual decrease In 
geold height of about 10 meters over a 1000 kilometer Interval. All of 
these tracks show both the deep low over the trench axis and the adja- 
cent high above the Island arc. Tracks G4263, G3681, and G3138 all have 
the gradual decrease behind the arc. Track G4806, however, has a slight 
Increase and then a step to a higher level. This step In geold height 
occurs at the continental margin southwest of Alaska. This same track 
also passes over the Emperor Seamounts and shows the characteristic high 
and flanking lows (Watts, 1979). Over the Kuril Trench track G3509 has 
the same features as the Aleutian Trench. 

An appropriate area to study deep-sea trench anomalies Is In the 
region of the Philippine Sea. This Is the location for the Mariana and 
Izu-Bonin, Philippine, Japan, and Ryukyu Trench systems. At the eastern 
boundary of the Philippine Plate the Mariana and Bonin Island arc-trench 
system has an Interesting sequence of geold undulations. Above the Izu- 
Bonin Trench tracks G4824 and G4861 In Figure 3 have an Increasing geold 
height as the trench axis 1s approached from the seaward side. Directly 
over the axis of the trench there Is a minimum of approximately 5 meter 
amplitude. This Is followed by an 18 to 25 meter increase to the high 
over the Bonin Islands. Adjacent to this Bonin high there Is another 
minimum which Is about 6 to 8 m In amplitude. Landward of the arc the 
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geold height remains fairly constant vrlth a gradual Increase. Further 
south above the Mariana Trench the Increase In geold height seaward 
of the trench Is more evident due to lengthier tracks of data with 65323 
and 64841. With these two profiles the Increase In geold height Is 
35 meters over 1500 ion distance. Over the Mariana Trench axis there Is 
a 10 meter mlnlimjm followed by a steep 20 iMter Increase landward of the 
axis. Above the Parece Vela Basin the geold Is fairly constant or has a 
gradual Increase. 

In the western Philippine Sea the plate boundary Is marked by the 
location of the Philippine Trench. 6E0S-3 geold profiles over this 
region are Illustrated In Figure 4. In these tracks, as they approach 
the trench axis from the West Philippine Basin, the geold gradually 
Increases In height by up to 20 m over 1800 km. Over the Philippine 
Trench Itself there Is a 5 to 15 m low, followed by a steep 15 to 25 m 
Increase to a geold high landward of the trench. Track 61694 passes 
over the Halmahera Ridge which has a 5 m high and two flanking lows. In 
the gravity field over this ridge there Is a +36 tngal anomaly and an 
adjacent -187 ixgal low (Watts. 1976). 

Some of the clearest examples of the geold ancxnaly at deep sea 
trenches occur In the region of the Tonga and Kermadec Trench system. 
These are Illustrated In Figure 5. Beginning In the Southwest Pacific 
Basin the geold Increases In height as the trench Is approached. This 
Is an Increase of almost 40 meters over a 2800 km distance. Similar to 
the other trench axes discussed previously, there Is a 10 m low over the 
axis of the Tonga- Kermadec Trench. This low Is followed by a steep 
Increase of up to 30 m to the high region over the South Fiji Basin. 

For 6E0S-3 track 62445 the geold Immediately begins to decrease landward 
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of the trench axis; from Figure 1 It Is seen that the track almost 
iRRiedlately Is back over the Pacific Basin. For tracks 61989 and G1477 
the geold decreases when the Pacific Basin Is similarly reached. Both 
tracks 62018 and 61590 cross over the South FIJI Basin and then have a 
local high above the FIJI Plateau. Track 61590 has a slight minimum 
above the South New Hebrides Trench. 

On the eastern side of the Pacific Ocean the Peru-Chlle Trench Is the 
location for subductlon of the Nasca Plate. This trench also has a 
gradually Increasing geold height as the trench axis Is approached from 
the seaward side. At this trench, however, the Increase Is very slight- 
only about 10 to 20 m over 1500 km (Figure 6). Again, as seen In the 
other regions, there Is a 5 to 7 m low over the trench axis followed by 
a 10 to 15 m step at the continental margin. Due to the presence of 
land, no altimeter measurements were made landward of the trench. 

Outside of the Pacific Ocean the only major trench system is In the 
Indian Ocean, the Java Trench. This region has the steepest geold 
anomalies seaward of the trench. At this location there Is a slope of 
about 55 m over a 1500 km distance (Figure 7). Above the trench Itself 
the low Is not as pronounced as In other regions. Here It Is only about 
2 to 5 m deep. Following the trench low there Is a dramatic 30 m 
Increase to the high region over the Java Sea. This Increase in geold 
height continues behind the arc at about the same rate as seaward of the 
trench . 

After observation of GEOS-3 altimeter data over these active 
subductlon zones, an Important conclusion Is possible. Over all major 
deep-sea trenches there Is a distinct pattern In the geold. Beginning 
anywhere from 1000 to 3(X)0 km seaward of the trench, the geold Increases 
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In height until Just prior to the trench axis. This very long Increase 
has a slope of 7 meters per 1000 km (1.4 seconds of arc) at minimum to 
37 meters per KXK) km (7.6 seconds of arc) at maxiimjm. In shape this 
Increase seaward of the trench Is In s«ne cases concave downwards 
(Aleutian Trench); concave upwards (Java Trench); or fairly linear 
(Tonga Trench). Directly over the trench axis there Is a narrow mlnlimjm 
In the geold. This low Is from 100 to 200 km In width and In depth 
anywhere between 3 m and 15 m (measured from the maxinum on the seaward 
boundary). From the base of the trench axis low there Is a steep 
Increase of up to 30 m to a high region landward of the tr: ch axis. In 
some cases there Is a distinct geold high above the Island arc. This Is 
most obvious with the Bonin Islands and also the Aleutian Arc. Landward 
of the Island arc the geold Is somewhat variable; behind the Aleutians 
there Is a gradual lowering of geold height. However, over the Java 
Sea the geold continues to rise. In other regions the geold remains 
fairly level landward of the trench. These features then are chaiacter- 
Istlc of the geold In the vicinity of deep-sea trenches; a grad v\1 long 
Increase as the trench Is approached from the ocean basin, a narrow low 
over the trench axis, and a relative high but with variable slope 
landward of the trench. 
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Gravitational Field In Vicinity of Deep Sea Trenches 

From the GEOS-3 altimeter measurements the characteristic geoid 
anomaly over deep sea trenches has been ascertained. In previous work 
at the ocean surface marine gravimeters have always been utilized to 
measure the gravitational field. These former studies have defined the 
free-air gravity anomaly which is characteristic of deep sea trenches. 

An early study of the Tonga Trench mapped the major features of the 
free-air gravity (Talwani, et al., 1961). In this work pendulum meas- 
urements were used to detect the large negative anomalies over the 
trench. Above both the Tonga and Kermadec Trench the free-air gravity 
has a minimum 100 to 200 km in width and -200 mgal in amplitude. Landward 
of the trench over the Tonga Ridge positive gravity anomalies, 100 to 
150 mgal in amplitude, were described by Talwani, et al. (1961). This 
region has a representative example of the typical free-air gravity 
anomaly above the trench and island arc. A gravity and bathymetry 
profile from Eltanin 29 is illustrated in Figure 5. For this profile 
there is a -220 mgal low over the Kermadec Trench and an adjacent 50 to 
100 mgal high above the Lau Ridge and Le Havre Trough. Other regions 
have similar gravity profiles, although at the Java Trench thore is a 
double gravity low with an intermediate high (Figure 7). 

Another common feature of the gravitational field is the existence 
of a positive free-air gravity anomaly seaward of deep sea trenches. 

This has been designated the Outer Gravity High by Watts and Talwani 
(1974). Originally observations of a bathymetric rise seaward of the 
Kuril Trench led Hanks (1971) to suggest horizontally compressive 
stresses occur in this region. According to his hypothesis, these 
stresses acting normal to the trench axis cause elastic deformation 
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of the lithosphere rrsultlng in an outer bathymetric rise. This concept 
was utilized by Watts and Talwani (1974) to explain the occurrence of the 
Outer Gravity High. In their model the elastic flexure of the litho- 
sphere causes the crust to bend upwards thus resulting in a positive gravity 
anomaly. Gravity and bathymetric profiles from Verna 20-06 in Figure 2 
show the outer topographic rise immediately seaward of the Aleutian trench 
and the associated positive free-air gravity anomaly 50 mgal in amplitude. 

A final feature of the gravitational field is the gravity anomaly 
landward of the trench. Above marginal basins the free-air values are 
near zero or slightly positive. As the crust in this area has a normal 
oceanic thickness yet occurs at shallow depths relative to the deep 
basin, it is possible to infer the existence of low density upper mantle 
(Karig, 1971). Beneath the Fiji Plateau Solomon and Biehler (1969) 
discovered a mass deficiency equivalent to a 1% density decrease to 
100 km depth located within the mantle. Yoshii (1972) calculated that 
low density mantle material occurs beneath the Sea of Japan. And behind 
the Izu-Bonin Trench, Segawa and Tomoda (1976) also calculated lower 
average mantle densities. What all of these studies indicate is the 
presence of lower mantle densities behind deep-sea trenches, when compared 
to the deep ocean basins. 

Combining all of these observations together, one can describe the 
characteristic free-air gravity in the vicinity of deep sea trenches. 
Immediately seaward of some trenches there is often an Outer Gravity 
High. This is about 50 mgal in amplitude and several hundred kilometers 
in width. Directly above the trench axis there is a narrow low, usually 
-150 to -250 mgal in amplitude. In the island arc region there are positive 
anomalies over the topographic highs with amplitudes of more than +200 
Landward of the island arc the free-air gravity anomaly is often close >, 
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zero or slightly positive. However, after elimination of crustal effects, 
lower average mantle densities are present when conpared to the deep 
ocean basin region. 
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Combination of Gravity and Geoid Measurements 

Technically, measurement of either the gravity field or geoid gives 
all the information possible about density variations within the 
Earth. Only one type of measurement is necessary. In a practical 
situation, due to the presence of noise, both altimetry and gravimetry 
measurements are useful for studying the density heterogeneities; they 
provide complementary information. To see why this is true, consider 
the nature of the geoid and gravity field. As observed in the previous 
figures (2 to 7), the geoid mostly has gradual undulations. Primarily 
it has long wavelength information. Free-air gravity anomalies, 
however, contain many fine details; in the gravity field short wave- 
length information is most evident. Due to the long wavelength char- 
acter of the geoid and the high wavenumber nature of the gravity field, 
for wavelengths longer than a certain number a radar altimeter will 
record the gravitational field with a higher signal to noise ratio; but 
at shorter wavelengths a marine gravimeter will perform better. This 
is illustrated in a quantitative manner by Chapman and Talwani (1979). 
When using data from both an altimeter and gravimeter the best possible 
information can be obtained at both shorter and longer wavelengths. 

For the Earth's gravitational field quite often long wavelength 
anomalies are ascribed to deeper structures while shorter wavelength 
features are explained by the more shallow structure. For example, 
gravity anomalies less than several hundred kilometers wide are often 
modeled by crustal structure as in the study by Talwani et a1. (1961). 
Long wavelength undulations typically are explained by phenomena within 
the mantle. Bott (1971) suggested the mantle transition zone was the 
cause of harmonics for wavelengths longer than degree and order 8. 
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This practice of relating longer wavelength anomalies Is explained 
by two principles. If a deep body has detailed structures and consequent 
high frequency gravity anomalies locally, these details will be lost due 
to attenuation with depth - only long wavelength Information will remain. 

A second principle Is that long wavelength gravity anomalies cannot be 
maintained by shallow structures. They would generate stress differen- 
ces In excess of the breaking strength of the shallow material. Com- 
bining these principles with the previous discussion. It Is possible to 
see that marine gravity probably gives the best constraints on the shallow 
densities while radar altimetry Is most useful for determining the 
deeper density heterogeneities. When data from these two sources are 
combined, there should be useful Information about both the deep and 
shallow lateral density variations within the Earth. 

In using both geold and gravity data to constrain Earth models, 
a problem arises with the gravity values. Very often in the region 
landward of the Island arc the free-air gravity anomaly is close to 
zero, yet the geold 1s high. This would appear to make interpretation 
quite difficult. It must be realized though that the free-air gravity 
anomaly Is caused by both mass inhomogeneities and changes In height 
of the geoid. Any effect due to geold height differences is often 
termed the 'Indirect Effect'. To Insure that gravity values are only 
caused by anomalous mass, a geold correction term should be added to 
the free-air gravity. These geoid corrections in marine gravity are 
discussed by Bomford (1962) and by Chapman and Bodine (1979). After 
correction for this factor both gravity ai;d geoid values may then be 
properly interpreted. 
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To give a concrete example of these Ideas, we utilize a previously 
published density model. In a thorough study of the Izu-Bonin Trench, 
Segawa and Tomoda (1976) used the free-air gravity anomaly and seismic 
refraction data to infer the crust and upper mantle density structure. 

Their model had an excellent match with the observed free-air gravity 
values. A reasonable ^^uestion to ask then is can this model be used to 
calculate geoid heights in agreement with the observed values. To 
answer this we utilized their model to compute the geoid anomaly. The 
result is illustrated in Figure 8. Above the calculated geoid height, 
the observed values from GEOS-3 are shown. The computed anomaly com- 
pares poorly with the observed data. This density model cannot explain 
the geoid anomaly over deep sea trenches. This poor agreement is due to 
two factors. As the calculated geoid values disagree with the measured 
data, there are obviously incorrect densities in their model. From our 
previous discussion it may be suggested that these incorrect density 
values are within the mantle. This is, however, not certain. Another 
important factor is the use by Segawa and Tomoda of the free-air anomaly. 
They neglected the 'Indirect Effect'. Landward of the arc, a value of at 
most 12 mgal should have been added to the free-air anomaly prior to inter- 
pretation. We have used this example to indicate that the use of geoid 
data in conjunction with gravity values provides a new and more powerful 
constraint on density models. The study of the Izu-Bonin Trench by 
Segawa and Tomoda (1976) was one of the most complete and excellent 
studies of a trench region; their model had very good agreement with the 
data they utilized. At this time, however, new data on geoid heights 
are available which would provide additional information to constrain 
their density model. 
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Possible Causes of the Geoid Undulation 


Due to the inherent ambiguity in gravity interpretation, it is 
always necessary to first consider the relevant physical processes for a 
given region. An understanding of these processes enables the range of 
possible density models to be greatly restricted. Existence of a posi- 
tive regional free-air gravity anomaly landward of active trenches has 
been previously known due to satellite based gravity models (Kaula, 

1969). A primary suggestion as to the cause of this anomaly is an 
excess mass due to descending lithosphere. While the slab penetrates 
into the asthenosphere it gradually heats up due to thermal conduction. 

As the rate of heat transfer is quite slow compared to the rate of sub- 
duction, the descending slab is in many cases cooler than the surroun- 
ding asthenosphere. Due to such thermal differences, the slab is denser 
than the adjacent material. This is an excess mass and causes a posi- 
tive gravity and geoid anomaly (Griggs, 1972; Hatherton, 1970 and 1969). 

In the calculations by Griggs the density contrast between the slab and 

3 

surrounding asthenosphere was at some depths greater than .2 gm/cm . 

This generated a gravity anomaly of +100 to +300 mgal in amplitude, 
depending upon the model parameters. Griggs also specified an undefined 
fonii of regional compensation. In a study of the crust and upper mantle 
beneath the Chile Trench, Grow and Bowin (1975) utilized density contrasts 
between .024 and .05 gm/om^ for the descending lithosphere. For the 
Izu-Bonin Trench, Segawa and Tomoda (1976) used a density contrast of 

3 

.05 p/cm for the slab. In a discussion on the gravity effect of 
downgoing slabs Watts and Taiwan! (1975) concluded the total o*',”'ity 
effect is limited in extent to the island arc-trench region and suggested 
that the regional positive anomaly is due to some other mass inhomogeneities. 
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In order to analyze the possible effect of a denser slab, a simple 
calculation has been made utilizing a two-dimensional technique (Chapman, 
1979). The hypothetical slab extends from 100 to 400 km depth where 

3 

there Is a density contrast of +.04 gm/cm with the surrounding mantle 
(Figure 9A). For this theoretical example there Is a gravity anomaly of 
+55 mgal and a geold undulation of 30 meters. The geold anomaly Is 
maximun over the center of the slab and decays downwards on either side. 
It should be noted that for all these two-dimensional calculations the 
absolute level of geold height Is arbitrary. 

Another possible suggestion as to the cause of the trench geold 
anomaly might be the outer bathymetric rise often found seaward of deep 
sea trenches. Such a model has been utilized to explain the free-air 
gravity seaward of the trench axis (Watts and Talwani, 1974). In the 
study of Watts and Talwani the bathymetric rise was considered to be the 
result of elastic flexure of the lithosphere due to regional stress 
fields. In Figure 9B such a possible model Is Illustrated. While the 
calculated gravity can explain the Outer Gravity High at trenches, the 
c(xnputed geold anomaly Is only 5 m high and a few hundred kilometers 
wide. As the observed geold anomaly at deep sea trenches Is much 
wider and larger In amplitude, the outer bathymetric rise would explain 
only a minute portion of the observed geold undulation. 

One plausible mechanism to generate the geold high Involves density 
changes within the mantle. Gravity studies In the past have Indicated 
a need for lower mantle densities (relative to the deep ocean basin) 
landward of the trench. This observation was made by Karig (1971) after 
studying seismic and gravity results at marginal basins. He noticed 
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POSSIBLE CAUSES OF GEOID UNDULATION 


DESCENDING LITHOSPHERIC SLAB 


outer bathymetric rise 



VARIABLE thickness LITHOSPHERE 
FULLY COMPENSATED 


VARIABLE THICKNESS LITHOSPHERE 
PARTIALLY UNCOMPENSATED 



^LANDWARD 


VARIABLE DENSITY LITHOSPHERE 
FULLY COMPENSATED 




VARIAK.E DENSITY LITHOSPHERE 

partially uncompensated 







Figure 9 






free*a1r values of gravity were near zero above these basins. Seismic 
results showed the crust In these areas to be of standard oceanic 
thickness, and yet to occur at shallower depths relative to the deep 
ocean basin. Arbitrarily assuming deep ocean basin to be normal, behind 
the trench there are higher densities at shallow depths. Relatively 
lower densities occur deeper. This Is a dipole type, layered mass 
distribution. Using the expressions of Ockendon d Turcotte 
(1977), this type mass distribution Is seen to cause a positive geold 
anomaly above the marginal basin. 

There are two possible ways to have lesser average mantle densities: 
either lateral changes In thickness or density of a body. Figure 9 C 
Illustrates the first possibility. In the region landward of the arc 
water depths are shallower, yet crustal thickness Is everywhere constant. 
At these shallow depths there are higher relative densities landward of 
the arc, when compared with the deeper ocean basin. Compensation 
occurs by having lower density (3.35 gm/cm ) mantle material at depth. 

If compensation Is not complete then the geold and gravity high are 
accentuated. This case Is also Illustrated In Figure 9 C. Figure 9 D 
Illustrates the second possible type mass distribution. Compensation 
occurs by having lower density material throughout the mantle to 100 km 
depth. An Important observation can be made by comparing Figures 9 C 
and 9 D. Both of these models cause virtually Identical geold and 
gravity values. Thus utilizing this type of data, the two types of 
possibilities are Indistinguishable. Additional constraints would 
be necessary to differentiate between the two causes. 

One final process which might affect geold undulations Is the 
presence of phase transformations within the descending lithosphere. 

Due to changes In pressure and temperature In the subducting slab, the 
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mantle material will undergo changes In phase at various depths (Ringwood, 
1976). Due to thermal differences, the transition at 400 Ion would 
migrate upwards In the slab with a consequent .2 gn/m Increase In 
density. If the slab Is considered a separate body, then It will 
conserve mass and the volume of the slab will decrease due to such phase 
changes. Decrease In the slab volume Insures that the total excess mass 
Is constant. After phase transformation the calculated geold and gravity 
anomaly Is virtually Identical to the values calculated prior to the 
change of phase. This Is because the total excess mass Is unchanged. 

Thus any occurence of phase transformations within the descending 11tho> 
sphere has no effect on the calculated geold and gravity anomaly. 

There Is of course no harm In Including phase transitions within the 
descending lithosphere, as many authors have done. 

After considering all of these physical processes, several obser- 
vations are possible. Any effect of the downgoing slab may be signi- 
ficant. However, the calculated anomaly Is maximtxn over the center of 
the slab and decays to either side. This calculated shape 1$ dissimilar 
to the observed geold undulation over deep sea trenches. The slab 
effect may be a contributory cause but cannot be the only mechanism 
which Is Involved. There must be some physical process to generate a 
generally high geold anomaly landward of the trench axis. A probable 
mechanism which will explain the geold high behind the arc Is a shallow 
crust either partially or totally compensated at depth by a mantle 
with variable thickness or density. Any presence of an outer bathy- 
metric rise will have only a minor effect on the calculated geold; phase 
transformations within the descending slab will have no effect on the 
calculated geold undulation. 
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Study of the Tonqa-Kermadec Trench 

Any attempt to ascertain the cause of the geold anomaly over deep- 
sea trenches Is best made by the explicit study of a given region. In 
the vicinity of the Tonga-Kermadec Trench the characteristic geold 
undulation Is quite evident and suitable data exists to study this 
area. One of the earliest studies of the Tonga Trench utilized pendulian 
gravity measurements aboard a submarine In combination with seismic 
refraction data. In this study by Talwani et a1. (1961) a crustal 
thickness of 36 km was calculated for the Tonga-Kermadec Ridge. As 
part of a comprehensive study of Melanesia extensive two-ship seismic 
refraction results were reported by Shor et a1. (1971). This latter 
work Is very useful because the seismic results give constraints on 
the thickness of the crust and sediments landward of the trench axis. 

In this present study gravity data were utilized which had been 
obtained on Eltanin cruise 29 from June to August 1967. Gravity 
measurements were made with a Graf Askania Sea Gravimeter Gss 2 
mounted on an Anschutz stable platform with an electrically erected gyro. 
Cross coupling errors were corrected by the use of an analog system. 
During this cruise a satellite navigation system In combination with 
the ship's electromagnetic log was used to provide accurate locations. 
Navigation, bathymetry, gravity and seismic profiles from Eltanin 29 
are presented In a report by Hayes et a1. (1972). 

Profiles of bathymetry, gravity and geold data across the Kermadec 
Trench are Illustrated In Figure 10. The gravity values are the 
gravitational attraction of anomalous mass. This is obtained by 
addition of a qeoid correction term (Chapman and Bodine, 1979) to 


161 



the free-air gravity anomaly. The effect of this correction is to 
increase the gravity values landward of the trench by 18 mgals maximum. 
Geoid measurements are from the GEOS-3 satellite. Both gravity and 
geoid values are referenced to a best fitting ellipsoid with flattening 
1/298.255. 

With the gravity data there is a 210 mgal high over Osborne Seamount 
(Part of Louisville Ridge) immediately seaward of the trench axis. 
Directly above the Kermadec Trench the gravity minimum is -210 mgal, 
while over Lau Ridge and Le Havre Trough gravity values average +70 
mgal with peaks up to +140 mgal. Landward of the arc system there is 
a regional high of about +35 mgal above the South Fiji Basin. With 
data from GEOS-3 track G1477 the geoid anomaly characteristic of deep- 
sea trenches is evident. Beginning 1600 kilometers seaward of the 
trench axis the geoid is +8 meters; gradually it Increases in height 
to +36 meters. Over the Kermadec Trench there is a low region in 
the geoid which is at a +27 meter level. Landward of the trench 
the geoid increases in height up to +56 meters, it remains at this 
high level above the South Fiji Basin. 

Any model of this region must be able to explain the general 
characteristics of the gravity and geoid. The density model must 
be able to produce a calculated geoid high landward of the arc, a 
local minimum above the trench axis and a gradual decrease in geoid 
height progressively seaward of the trench axis. From the bathymetry 
it is evident that the South Fiji Basin, Lau Ridge and Le Havre Trough 
are shallower than the Southwest Pacific Basin. These shallower 
depths must be at least partially compensated otherwise there would be 
much higher gravity values landward of the trench. Of primary interest 
then is where the lower average densities must occur. 
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A standard assumption in gravity modeling is that variations in 
crustal thickness can be the mechanism of compensation. In our model 
this was not the sole mechanism. Our final crustal model is illustrated 
in Figure 11. The gravity effect of this is also shown in comparison 
with observed data. Quite obviously this model produces gravity values 
much greater than observed data for the region landward of the trench 
axis. This calculation indicates the necessity for lower average 
densities to be present in the mantle landward of the trench. A similar 
conclusion has been reached by other authors for different trench regions; 
for example, Segawa and Tomoda (1976) in the Izu-Bonin region. 

This model of the crustal structure was obtained by consideration 
of seismic refraction data from Shor et al. (1971) to help constrain 
crustal thickness. As their measurements were obtained about 60 nm 
to 300 nm south of the Eltanin 29 track, they were projected along strike 
to the ship profile. For the Southwest Pacific Basin no refraction data 
existed and a standard oceanic crust was assumed. In the standard 

3 

model there was 5.5 km water, 2.0 km upper crust with density 2.7 gm/cm , 

3 

4.0 km lower crust with a 2.9 gm/cm density overlaying mantle material 

3 

with density 3.4 gm/cm . In the Southwest Pacific 'sin any deviations 
from the standard model were fully compensated in an Airy manner. After 
inclusion of mantle densities, crustal structure landward of the trench 
was slightly adjusted until calculated gravity values agreed with 
observed data. Even though the refraction data was not perfectly 
coincident with the ship track, our final crustal model is in good 
agreement with seismic results. Table 2 has a comparison of the density 
model and seismic refraction values of Shor et al. (1971). 

For both this crustal model and the entire model including the deeper 
mantle a spherical Earth method was used. Johnson and L’tehiser (1972) 
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Figure 11 
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developed a technique for calculation of gravitational potential and 
attraction due to a spherical cap with polygonal sides. Geoid height is 
computed by dividing anomalous potential by normal gravity. In our 
utilization of this method the body is a input in cross section. Each 
segment is a spherical lamina which extends perpendicular to the cross 
section a certain distance. In our model the width was 2000 kilometers. 
Modeling was initially done assuming a plane two dimensional earth, 
later modeling or>ly used spherical earth calculations. It was discovered 
that due to neglect of curvature the plane earth models gave erroneous 
values seaward of the trench. To insure that anomalous mass was being 
used, the gravitational effect of a standard Earth was subtracted. The 
standard model contained standard oceanic crust as previously defined and 

3 

uppermost mantle 88.5 km thick with density 3.40 gm/cni . Only hori- 
zontal differences in density are important thus any effect of adiabatic 
compression was neglected. Our model thus contains no information on 
vertical density contrasts. All that is important are the lateral 
density differences. 

Inclusion of mantle structure completes our final density model as 
shown in Figure 12. Crustal structure is identical to that which is 
illustrated in Figure 11. This model contains structures which were 
defined by seismic research. Seismicity studies have shown the presence 
of the Benioff zone (Sykes et al., 1969), while seismic velocity and 
attenuation anomalies (Mitronovas and Isacks, 1971; Oliver and Isacks, 
1967) indicate this zone is caused by descending lithosphere. In our 
model * location of the descending lithospheric slab is based on the 
most recent seismicity study by Billington (1979). Thickness of the 
slab was constrained at 55 km down to a depth of 250 km and at 40 km for 
depths to 600 km. Studies of attenuation of shear waves indicate the 
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Figure 12 



presence of a low Q zone above the slab near Tonga (Oliver and Isacks, 

1967). Possibly this Indicates an anomalous zone, thus a region of low 
density material was Included between 65 km and 180 km depth beneath the 
Island arc. This anomalous zone may extend to shallower depth or also 
further landward. Because there was no definite seismic evidence concerning 
Its upper limit, this was not done. 

Our procedure In modeling was to first use seismic refraction data 
to constrain crustal structure. For the mantle location of a descending 
slab and low Q Zone were then fixed, constrained as well as possible by 

seismic studies. Density of the uppermost mantle was assumed to be 

3 3 

3.4 gm/cm while underlaying mantle was fixed at 3.35 gm/cm . For the 

mantle region densities within the slab and low Q Zone were varied, as 

was the thickness of the uppermost mantle landward of the slab. Minor 

adjustments In crustal structure were also made until calculated values 

of geold height and gravity agreed with the measurements. 

Mantle structure In this final nrodel has the following charac- 
teristics. Depth to the lower density mantle landward of the trench Is 
35 km less than beneath the Southwest Pacific Basin. Thinning of the 
uppermost mantle body Is the mechanism for generating a mass deficit 
within the mantle. This Is because density of the lower density mantle 

3 

Is .05 gm/cm less than density of the uppermost mantle body at similar 
depth. These lower average densities are the compensating mechanism for 
the higher average densities at shallow depths landward of the trench. 
However, this model Is not fully compensated. At a depth of 100 km 
there Is a difference In pressure of 60 bars beneath the South Fiji 
Basin compared to the Southwest Pacific Basin. There Is an excess of 
mass amounting to 7.3 x 10^ gin/cm^ In a column beneath the South Fiji 
Basin relative to the Pacific Basin. This mass excess would correspond 
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to a +31 mgal anomaly over the South Fiji Basin, assuning a Bouguer 
approximation. 

In this model the effect of the slab Is relatively small. Compared 

3 

to the surrounding mantle the slab has a density contrast of +.02 gm/cm . 
During our modeling procedure It was discovered that If the density 
contrast of the slab was +.05 gm/cm or larger there would be no calculated 
low In the geold over the trench axis, although there would still be a 
-190 mgal m1n1mimi above the axis. Use of the geold data constrains 

3 

the density contrast of the slab to be less than .05 gm/cm relative to 
the surrounding material. In order to match the observed data, however, 

3 

the density contrast was only +.02 gm/cm . Gravitational effect of the 
slab Is Illustrated In Figure 13. Both geold and gravity are maximum 
above the slab with values of +13 m and +24 mgal respectively. Presence 
of lesser densities In the low Q zone negate the effect of the slab. Due 
to both bodies the total effect Is -8 m and -38 mgal at maxlmun. As 
regards the effect of the slab, several conclusions are possible. There 

3 

is a maximum density contrast of +.05 gm/cm for the slab, although In 

3 

our model It Is required to be a lower value, +.02 gm/cm . Calculated 
geold heights and gravity values due to the slab explain only a small 
portion of the observed anomaly. Presence of lesser densities In the 
low Q zone counteract the Influence of the slab and result In a total 
effect which Is slightly negative. As the net effect of both the slab 
and low Q zone explains only a small portion of the observed anomaly, It 
might be possible to construct a total model excluding those bodies. 

Such a possibility was not tested however, because there Is definite 
seismic evidence for these bodies. But It Is clear that from gravity 
and geold evidence alone, the evidence of a dense slab cannot be Inferred. 
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GRAVITATIONAL EFFEC^OF DESCENDING SLAB 
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Suwnarv and Conclusions 

After observing GEOS-3 radar altimeter data, a characteristic geold 
anomaly has been Identified which occurs over all major deep-sea trenches. 
Seaward of the trench there Is a gradual Increase In geold height with 
higher values near the trench and lower geold heights progressively seaward. 
Slope of this very long Increase ranges between 7 meters and 37 meters per 
1000 km of distance. This gradual slope begins anywhere between 1000 and 
3000 km seaward of the trench. In some regions the slope Is concave downward 
as at the Aleutian Trench. Other areas have a slope which Is fairly linear 
(Tonga Trench) or concave upward (Java Trench). Directly above the axis of 
trenches the geold forms a narrow low from 100 to 200 km 1n width and 
anywhere between 3 m and 15 m In depth. Beginning at the base of the trench 
axis low there Is a steep Increase In geold height of up to 30 m to a high 
region landward of the trench axis. Over some Island arcs, for example the 
Bonin Islands or the Aleutian Arc, there Is a distinct r^'^old high. This Is 
between 5 and 10 m In amplitude. Landward of the trench axis the geold has 
a generally high level relative to the deep sea basins. In this region 
however, slope of the geold Is variable. Behind the Aleutians the geold 
has a slow decrease progressively landward, while over the Java Sea there 
Is a progressive Increase. In other regions the geold remains fairly 
constant landward of the trench. 

In utilizing data to study the Earth It Is best to combine both geold 
and gravity measurements. Gradual undulations of the geold constrain 
densities deep within the Earth while gravity anomalies help determine 
crustal density. In using both types of data it is important to correct the 
free-air gravity anomaly for undulations of the geold. Addition of a geold 
correction term insures that gravity values are equal to the attraction of 
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the anomalous mass. 

A density model of the Tonga«Kermadec region has been constructed 
which explains both the observed gravity and geold values. Occurence 
of both a geold and gravity high landward of the trench Is caused by a 
dipolar type mass distribution. Landward of the trench axis, water 
depths are shallower when compared to the deep ocean basin. Thus at 
shallow depths there are higher densities relative to the Southwest 
Pacific Basin. An abnormally thick crust Is not the mechanism for 
compensating these higher relative densities. Rather, In our model, 
compensation occurs by thinning of the uppermost mantle body. Beneath the 
South Fiji Basin lower density mantle material occurs at a depth 35 km 
shallower than beneath the Southwest Pacific Basin. These lower relative 
densities do not provide full compensation. At 100 km depth, pressure 
beneath the South Fiji Basin 1s 60 bars greater than beneath the Southwest 
Pacific Basin. 

Perhaps one could associate the uppermost mantle body (density equal 

3 

to 3.4 gm/cnr) with lithosphere and the less dense, underlying mantle 
body, with a^thenosphere. Using this terminology, partial compensation 
occurs within the mantle by a thinning of the lithosphere beneath the 
South Fiji Basin. According to Watts et a1. (1977) age of the ocean 
floor In this region Is approximately 30 m.y.b.p. If this were normal 
oceanic lithosphere, depth to the asthenosphere would be expected to be 
approximately 65 km (Forsyth, 1977). This coincides with our model and 
Is an Interesting observation, but most likely fortuitous. 

In explaining the geold anomaly any effect of the descending 

lithospheric slab is limited. Density contrast for the slab cannot 
3 

exceed .05 gm/cm , but in our model is only .02 gm/an . In our model 
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gravity values due to the slab are mgal at maximum. 

In the characteristic geold anomaly at deep sea trenches there Is 
an Increase In geold height as the trench Is approached from the seaward 
side. This gradual sloping of the geold seaward of the trench Is 
entirely explained by the mass distribution landward of the trench axis. 

No anomalous mass (relative to landward regions) beneath the Pacific Basin 
Is required to explain the gradual slope In the geold. 

During our modeling. It was found that It was quite difficult to 
match both gravity and geold values. This Is a qualitative observation 
that they provide complementary Information. Also the geold was quite 
sensitive to changes In mantle structure. If the depth to the lower 
density mantle material landward of the arc was altered by 5 km, there was a 

10 m change In calculated geold height. Yet In gravity this caused only 

a 10 mgal difference. 

As In any study of the gravitational field there are limitations to 
our model. For exan^le landward of the trench we cannot distinguish 
between thinning of an uppermost mantle body or lesser densities within 

the uppermost mantle. A thinner mantle body Is the mechanism we utilized 

to produce a mass deficit at depth, but a seismic study would be nece- 
ssary to confirm this conclusion. Similarly for the descending slab a 
qualification Is necessary. For a slab thickness of 55 km thinning to 

3 

40 km. It was found that the density contrast could not be .05 gm/cm or 
greater. Otherwise the geold low above the trench axis could not be 
computed. If however a thicker slab were utilized, the density contrast 
would decrease. This data would not be able to distinguish between 
different thickness of descending lithospheric slabs. Lastly, In our 
model density within the mantle was assinned to have no lateral variation. 
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This Is an untested assumption, but perhaps seismic surface wave studies 
In the ocean basins could detect any such variability within the 
mantle. 

Study of the Tonga-Kermadec region raises Interesting questlems. 

Can the sloping geold anomaly seaward of all trenches be explained by 
mass distribution landward of thv trench axis? As this anomaly Is 
variable In shape, perhaps mass heterogeneities are also required 
beneath the deep ocean basins. 

Another enigma Is the variability In geold height landward of the 
trench axis. Can a model be constructed which would produce a variable 
geold slope In this region? Our modeling seemed to suggest this would be 
possible by varying the horizontal extent of the dipolar type mass 
distribution. 

Another Intriguing lesult In our model Is the excess pressure at 
100 km depth landward of the trench. What are the forces required to 
sustain this model in equilibrium? Presumably any dynamic model which 
can generate such forces, would give rise to a geold and gravity ano- 
maly. In cur model any variation In density within the asthenosphere 
was neglected, yet such differences would be required to maintain our 
model In equilibrium. Sleep (1975) has suggested a model for flow 
beneath Island arcs. In his work topography was explained by stresses 
amounting to several hundred bars. It would be most interesting to 
construct a geophysical model capable of explaining any possible dynamic 
flow along with observed geold and gravity anomalies. 
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APPENDIX 


A) Finite Vertical Sheet 


N 


Z2 u ^2 

6 0 j logg^(x-x') + [(x-x')2+ (y-y‘)^ + (z-z')jj dz 

21 X] 


This Integral Is of the form 
logg fa+ [b^ + x^J^x 

9 9 0 

letting u ■ B‘ + x‘, restricting u. > o and integrating by parts we obtain 

■ (u^ - B^^log. (a + u) - Z' Cu^ - B^J ^du 

w/ a + u 

and after Integrating again 

. 0,2 - (a ♦ U) - C.2 - b2]’* 

6 r T 

+ a log« /Cu^ - B?I+ ui +Cb^ - a2)sin"^ -1 (b^ + au)| 

62 -a 2 ]’i I J 


or finally 

y'logg {a + + x gjdx = X logg ^a + 0,2 * X 

+ a log^ {x + fc2 + x^J+ ^ 

[b2 - a2feln-’ r.l/b2*a ^2 ^x2^ 
[la* [b2 + x3'^B| 
restricting x > o 

B) Two dimensional rectangle 

I2 


N 




(y - y') log- C(y - y'r + (z - z')^] -2y 


+ 2 I z - z 'J Tan 


-1 


(y - y') 

13 

1 

M 

cr 



yz 


^1 


dz 


restricting (z - z') > o (this is the same integral as used for the two dimensional 
sheet) and integrating 
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C) Two dimensional polygon 

For the ^th side of the polygon we have 

"l “ C(y-y')^ + (z-z')^] - 2 (y-y') 

+ 2 (z-z') Tan"^ |v-y ' | j dz 

and we Integrate along a linear line segment defined by z = m-f y + b-j ; 
assuming m^ and 1/m-j are not zero we obtain 
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L u-z') J 

where m-j ■ Zj - 


y\ - yi -1 

Bi » z^ - mi yi 

First part of the Integral is of the form log^ (a+bx+cx^) dx and is a 
standard integral. Third part of the integral is of the form y'x Tan"^ 
[(Cx + D)/x] dx. Letting u « (Cx + D)/x and integrating by parts, this 
becomes 



After the final integration we find for the third part 
^x Tan" [(Cx+D)/x] » 7 ^x^ Tan"^ [(Cx+0)/x] 

- C logg r (1 -h:^) x^ 2CDx + dM 

( 1 ^ 2)2 L .2 J 


+ D_x__ + XI^cVd^ Tan"^ [(Cx+D)/x]^ 
(1+C2) (1+C^)^ 




Thus formula 23 Can be integrated to obtain formula (24). 

In evaluating the line integral the following special cases are also 
important. 

Case (a). If m. = 0, then dz is zero and consequently ni is also zero. 


183 


Case (b). If 1/m^ « 0 then y Is constant along the line segment and the 
line Integra. (23) is written as: 

n^ • y* [(yry') ^09e C(yt-y')^ + (z-z')^] 

^1-1 

-2 (yi-y') + 2 ( 2 - 2 ') Tan"^ (yi-y') *^ dz 

(z- 2 ’) ' 

This Is Identical In form to the Integral for the two dimensional rectangle 
and Is solved In an Identical fashion. 

Case (c). If the extension of the line segment passes through the 
observation point (when B,j/m^ + y' ■ 0) then the line Integral (23) Is 
written as 


*•1 

n^ » logg [(l+l/m^^) 2 ^] - 2 2 /m^ 

+ 2 Z Tan”^ (1/mi)} 

assuming Z' * 0, Z > 0 

This Is easily Integrated to become equation (27). 

D) Three dlmenslc.nal polyhedron. For the j'-th line segment of the 

rth facet we have, after Integrating equation 36 by parts and rearranging 

ni+i 

Lj * Gp0i C n ('"jn+Cj) + C(">jn+Cp^ * 

9 .. V 


^1+1 
/' 2 
/’ ^ dn 

J n^+ 0 ^ 


'’l+l 

+ f Jli, V 
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'^ 1+1 





(-J»^ * d,1 



The first Integral Is Independent of the path of Integration, thus when 
Is surrmed over all line segments for a closed path this Integral will sum 
to zero. Consequently, we Ignore this term. The second Integral Is In a 
standard form. To solve the third Integral the substitution Is made. 

n ■ «1 ■ '3^ * 

(mj Cj) (x+iy 

This will transform the last Integral Into 

['"j»f*'3!l'^"^i'3 f .,,‘4 

which Is a standard form and can be Integrated. This same integral (36) has 
also been solved for gravity problems by Barnett (1976) and Paul (1974). 
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